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has  been  used  to  obtain  traction  measurements  for  a  MIL- L- 7808  oil.  A  full  matrix  of 
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thermal  theory.  Capacitance  measurements  obtained  at  higher  loads  and  speeds  indicate 


more  significant  load  variations  than  existing  theory  would  predict. 


The  partial  elastohydrodynamic  regime  has  been  studied  analytically  and  experimentally. 
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measured  electrical  resistance  and  predictions  obtained  with  the  asperity  interactions 
analysis  put  forth  in  Task  2  have  been  obtained. 
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SECTION  I 


INTRODUCTION 


At  the  conclusion  of  Phase  I  of  this  research  program,  a  compendium  of  the  state 
of  the  art  of  elastohydrodynamic  technology  was  compiled  and  a  computer  program 
was  prepared  to  predict  elastohydrodynamic  performance  characteristics  on  the 
basis  of  the  best  available  theory  and  experimental  data.  It  was  found  that 
considerable  data  were  lacking  in  areas  such  as  elastohydrodynamic  tractions, 
measurement  of  elastohydrodynamic  film  thicknesses,  and  general  behavior  of 
partial  elastohydrodynamic  performance. 

* 

A  rolling  disc  elastohydrodynamic  test  apparatus  was  built  and  instrumented  (1) 
as  a  result  of  the  findings  in  Reference  (2),  for  measurement  of  elastohydrodynamic 
performance  characteristics  including  tractions  and  film  thicknesses.  A  full 
matrix  of  traction  data  was  obtained  for  a  SP^E  polyphenyl  ether  during  Phase 
II  of  this  research  program  along  with  analytical  predictions  to  describe  the 
data.  An  asperity  interaction  computer  program  was  also  published  in  Refer¬ 
ence  1  as  a  vehicle  for  studying  the  partial  elastohydrodynamic  regime. 

During  Phase  III  of  this  program,  a  full  matrix  of  traction  data  was  obtained 
with  a  MIL-L-7808  oil.  Comparisons  made  between  this  oil  and  polyphenyl  ether 
showed  radical  differences  in  both  quantitative  and  qualitative  traction  be¬ 
havior  even  under  conditions  where  conventional  properties  are  relatively 
similar. 

Attempts  were  made  to  obtain  film  thickness  data  using  three  different  tech¬ 
niques;  namely,  optical,  capacitance,  and  X-ray  techniques.  Optical  data 
was  obtained  which  showed  good  detail  of  the  film  thickness  profile  at  lower 
loads.  Thes  ■  were  compared  with  capacitance  data  taken  with  steel  discs  under 
conditions  of  higher  loads  and  speeds.  Calibrating  techniques  were  developed 
for  relating  capacitance  measurements  to  film  thickness  ar>H  for  checking  out 
optical  measurements  under  elliptical  contact.  Full  dry  contact  elastic  Us- 
tortion  analyses  for  Hertz  contact  were  considered.  The  optical  data  taken  at 
lower  loads  indicated  relatively  small  load  sensitivity  as  is  predicted  by 
elastohydrodynamic  theory.  Capacitance  data  taken  at  higher  loads  in¬ 
dicated  much  stronger  load  effects  which  are  in  qualitative  agreement  with 
X-ray  data  given  in  Reference  3. 

Electrical  continuity  and  resistance  measurements  were  obtained  with  the  use  of 
an  asperity  contact  measuring  device  developed  by  MTI  and  referred  to  as  an 
"Asperitac."  This  device  provides  information  regarding  the  percentage  of  time 
the  electrical  resistance  falls  below  a  preset  discrimination  level  and,  hence, 
is  a  measure  of  the  relacive  extent  of  contact.  It  was  found  for  the  cases 
studied  that  no  significant  contact  occurred  in  the  higher  speed  regions  where 


Numbers  in  parentheses  refer  to  references  listed  at  the  end  of  this  report. 
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capacitance  and  traction  measurements  were  obtained,  and  in  order  to  get  incip¬ 
ient  contact,  speeds  well  below  the  100  in/sec  range  had  to  be  used.  In  general, 
based  upon  elastohydrodynamic  film  thickness  predictions,  this  would  be 
the  range  in  which  the  film  thickness  to  surface  roughness  ratios  are  on  the 
order  of  1.  Comparisons  were  performed  with  electrical  resistance  data 
and  resistance  predictions  based  upon  the  asperity  interaction  model;  however, 
conclusions  are  not  definite  as  a  result  of  uncertainties  in  the  distribution 
of  asperities. 

The  following  sections  deal  first  with  measurements  and  correlation  of  trac¬ 
tions  for  MIL-L-7808  fluid.  This  section  is  followed  by  three  sections  dealing  1 

with  film  thickness  measurements  obtained  by  capacitance,  optical,  and  X-ray  J 

techniques.  These  are  followed  by  the  results  of  the  asperity  interaction  1 

studies.  Appendices  are  presented  to  give  the  full  matrix  of  traction  data  \ 

1  along  with  the  various  lubricant  data  and  mathematical  analyses  that  have  been 

;  performed  in  relationship  to  the  traction  and  film  thickness  data.  \ 
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SECTION  II 

TRACTION  DATA  AND  ANALYSIS 


This  section  presents  a  summary  of  the  traction  data  collected  from  tests 
performed  on  a  MIL-L-7808  lubricant  along  with  a  ouasi-empirical  model  for 
predicting  tractions.  The  test  lubricant  physical  properties  are  included  in 
Appendix  I  for  reference  purposes.  The  test  rig  used  to  generate  the  traction 
data  discussed  here  was  reviewed  in  Reference  1. 

1 .  TEST  DATA 

The  test  rig  design  is  of  the  two-disc  variety  which  permits  the  determination 
of  traction  and  film  thickness  properties  of  most  lubricants  in  a  single  roll¬ 
ing  and/or  sliding  contact.  In  addition  to  accurate  control  and  readout  of 
rolling  and  sliding  speeds,  various  levels  of  temperature  can  be  maintained  on 
the  disc  machine  for  investigating  the  dependence  of  traction  on  lubricant 
temperature.  A  sliding  carriage  supporting  the  upper  disc  specimen  also  pro¬ 
vides  a  way  to  vary  the  load  in  order  to  study  loading  effects  on  a  lubricated 
contact. 

The  basic  geometry  of  the  test  disc  specimens  was  identical  to  that  used  during 
the  past  year's  work  completed  with  polyphenyl  ether.  Tlie  traction  data  was 
taken  from  two  crowned  disc  specimens,  each  with  a  three-inch  diameter  and  a 
thirty-six-inch  crown  radius.  Both  were  hardened  and  polished  AISI  52100  bear¬ 
ing  steel.  The  tested  hardness  of  each  was  between  fifty-eight  and  sixty 
Rockwell-C.  The  finish  of  the  polished  surfaces  of  the  rollers  was  approxi¬ 
mately  two  microinch  rms.  A  topographical  recording  of  five  parallel  tracings 
of  a  typical  roller  surface  used  for  the  tests  is  shown  in  Figure  1.  Each 
tracing  was  recorded  along  a  line  perpendicular  to  the  direction  of  rolling  on 
the  specimen  surface,  A  complete  discussion  of  test  surface  topology  as  rela¬ 
ted  to  asperity  interactions  of  the  test  specimens  is  presented  in  Section  VI. 

Full  traction-versus-slip  curves  for  both  positive  and  negative  slips  are  a 
direct  output  of  the  rig  instrumentation.  A  sample  of  four  typical  traction 
curves  read  out  directly  from  the  disc  machine  is  shown  in  Figure  2.  Each  of 
the  four  sample  curves  was  obtained  at  different  loads  but  constant  rolling 
speed  and  lubricant  temperature. 

Traction-versus-slip  curves  for  a  wide  range  of  contact  loading,  rolling  speed, 
and  lubricant  temperature  have  been  investigated  with  a  MIL-L-7808  lubricant. 
The  range  of  conditions  for  presently  available  traction  data,  including  the 
previous  year's  work  performed  on  polyphenyl  ether,  is  shown  in  Figure  3.  In 
order  to  cover  the  range  of  variables,  a  large  number  of  complete  slip  curves 
have  been  recorded  (to  date,  135).  Enough  combinations  of  loading,  rolling 
speed,  and  lubricant  temperatures  have  been  recorded  to  reasonably  describe 
the  traction  behavior  of  the  two  lubricants  tested. 


Since  a  large  number  of  complete  traction  curves  can  be  obtained  in  a  short 
period  of  time,  the  data  handling  and  storage  has  been  tied  into  a  computer 
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Fig.  1  Five  Parallel  Talysurf  Tracings  of  Roller  Specimen  Surface  Finish 


Fig.  2  Traction  Data  from  Rolling  Disc  Rig  for  Four  Loads 
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Fig.  3  Range  of  Presently  Available  Traction  Data  From  MTI  Computer 
Data  Storage 
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via  time-sharing  terminal.  The  curves  are  digitized  and  transferred  in  total 
to  the  computer  filing  system  at  the  present  time  by  hand.  The  two-fold 
symmetry  axes  for  each  complete  traction  curve  is  found  analytically  and  a 
single  positive  slip  curve  is  stored  in  a  computerized  traction  data  bank. 


Retrieval  and  analysis  of  the  data  is  highly  automated  since  the  present  time¬ 
sharing  facility  at  MTI  is  linked  to  an  onsight  plotter  which  can  display  any 
stored  curve  within  seconds.  A  sample  set  of  the  computer-plotted  curves  is 
shown  in  Figure  4,  which  shows  a  comparison  of  the  four  loads  previously  dis¬ 
played  as  raw  data  in  Figure  2.  The  scale  of  the  absicca  and  ordinate  can  be 
plotted  in  any  convenient  set  of  engineering  units  from  the  stored  files  of 
traction  data.  The  present  data  handling  system  allows  rapid  comparison  of 
gathered  data  taken  under  any  set  of  test  facility  conditions. 


The  computer-plotted  comparison  of  a  5P4E  polyphenyl  ether  and  a  MIL-L-7808 
oil  is  shown  in  Figure  5.  The  conditions  of  disc  geometry,  load,  and  rolling 
speed  are  identical.  The  lubricant  temperature  conditions  were  widely  separa¬ 
ted  in  order  to  match  the  viscosities  of  the  fluids.  Both  fluids  have  a 
pressure  coefficient  of  approximately  0.75  x  10“^  in^/lb.  One  curve  shows  the 
traction-slip  rate  for  polyphenyl  ether  at  260°F,  whereas  the  other  curve 
shows  traction  data  obtained  from  the  MIL-L-7808  oil  at  120°F.  Both  viscosi¬ 
ties  are  approximately  nine  centipoise. 


I 


* 

- 

5 


A  marked  difference  is  displayed  even  though  the  viscosity  and  viscosity-pressure 
coefficient  characteristics  are  similar.  Part  of  the  displayed  difference  can 
be  attributed  to  the  fact  that  the  temperature  coefficient  of  viscosity  for  the 
MIL-L-7808  oil  is  about  80  percent  greater  than  that  for  polyphenyl  ether.  How¬ 
ever,  this  is  insufficient  for  explaining  the  wide  discrepancy  revealed  from  the 
data. 

Another  Indication  of  the  fact  that  conventional  viscosity  characteristics  are 
insufficient  for  describing  elastohydrodynamic  tractions  results  from  the 
polyphenyl  ether  traction  data  obtained  by  MTI  which  indicates  tractions  for 
polyphenyl  ether  to  be  nearly  insensitive  to  temperature.  Conventional  viscos¬ 
ity  characteristics,  such  as  the  base  viscosity  and  viscosity-pressure  and 
temperature  coefficients,  vary  markedly  with  temperature.  The  lack  of  sensi¬ 
tivity  to  temperature  indicated  by  this  data  may  result  from  additives  intro¬ 
duced  into  the  polyphenyl  ether  (Battelle  data  indicates  considerable  tempera¬ 
ture  sensitivity  in  tractions  for  polyphenyl  ether  without  additives).  The 
effect  of  these  additives  does  not  show  up  in  conventional  viscosity-temperature 
data. 

The  present  test  lubricar. :  (MIL-L-7808),  however,  showed  a  strong  effect  with 
different  operating  temperatures.  Shown  in  Figure  6  is  a  display  of  the  temp¬ 
erature  effect  exhibited  for  one  constant  load  and  fixed  rolling  speed.  As  the 
temperature  increases,  the  maximum  transmittable  peak  traction  is  reduced. 

Figure  7  ahows  a  25  percent  reduction  of  the  peak  traction  was  observed  for  an 
85 °F  variation  in  temperature  at  the  speed  and  load  indicated.  One  additional 
set  of  traction  curves  taken  at  higher  rolling  speeds  is  plotted  in  Figure  8. 
Included  in  Appendix  II  for  reference  are  computer  plots  of  all  the  MIL-L-7808 
traction  curves  stored  to  date  for  retrieval  from  the  MTI  computer  data  bank. 
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Fig.  4  Typical  Load  Sensitivity  of  MIL- L- 7808  Traction  Versus  Slip  Data 


•8- 


HTl-13432 


TRRCTI0N  <L&> 


TRRCTIHN  ZDRTR 


BLIP  RRTE  <IP5> 

Fig.  5  Comparison  of  Disc  Traction-Slip  Data  for  Two  Different  Lubricants 
with  Similar  Viscosities 
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Peak  tractive  forces  varied  widely  as  a  function  of  the  applied  Hertzian 
pressure,  rolling  speed,  and  temperature.  Since  a  large  number  (over  200) 
of  complete  traction-versus-slip  curves  were  collected,  cross  plots  of  roll¬ 
ing  speed  versus  Hertzian  pressure  for  fixed  peak  tractive  forces  can  be 
made,  A  single  line  for  polyphenyl  ether  is  shown  in  Figure  9  for  a  peak 
tractive  force  of  four  pounds,  A  range  of  peak  tractions  greater  than  one 
pound  and  less  than  four  is  shown  for  the  MIL-L-7808  tested  in  Figure  9, 

2.  ANALYSIS  OF  TRACTION  DATA 

An  attempt  has  been  made  to  analyze  the  traction  data  and  correlate  it  with 
the  use  of  available  viscosity  data-  A  quasi-empirical  model  has  been  devel¬ 
oped  for  predicting  tractions  based  upon  hypothetical  but  reasonable  viscosity 
data. 

The  predicted  tractions  along  with  the  measured  ones  are  graphically  displayed 
in  Figures  10  through  15.  The  extremes  in  rolling  speed  and  lubricant  test 
temperatures  have  been  plotted  for  comparison.  For  the  range  of  test  condi¬ 
tions,  it  can  be  seen  that  satisfactory  predictions  have  been  obtained 
which  may  be  used  to  describe  MIL-L-7808  tractions  over  the  full  test  range  of 
loads,  rolling  speeds,  slip  rales,  and  temperatures. 

The  methods  that  have  been  used  for  obtaining  the  traction  correlation  are 
described  in  the  following  sections, 

a.  Film  Thickness  Predictions 

Any  hydrodynamic  traction  prediction  requires  a  knowledge  of  shear 
rate  as  input.  In  order  to  compute  the  shear  rate  from  the  known 
slip  velocity,  an  estimate  of  the  film  thickness  is  required-  This 
estimate  is  obtained  with  the  use  of  the  elastohydrodynamic  perform¬ 
ance  code  prepared  for  the  Air  Force  and  listed  in  Reference  2 , 

The  hydrodynamic  film  thickness  inputs  to  the  film  thickness  portion 
of  this  computer  program  are  the  thermal  conductivity  of  the  lubri¬ 
cant  and  the  viscosity  coefficients  a,  and  3  corresponding  to 
the  formula. 


e  p  +  (2  +  vp ) ( T  +  A60 


_i _ ) 

T  +  460' 
o 


where  T0  denotes  the  reference  temperature  in  degrees  Fahrenheit  taken 
here  to  be  the  lubricant  supply  temperature  and  ~0  is  the  viscosity  at 
temperature  T0  Values  of  ~0  and  2  were  determined  from  atmospheric 
viscosity-temperature  measurements  obtained  at  MTI  and  shown  in  figure 

16, 

Values  of  -  are  obtained  from  the  local  slopes  of  a  semi-log  plot  c: 
h  versus  1/TA  where  is  the  absolute  temperature  in  accordance  with 
the  relationship 


Peak  Traction  uepenueuee  on  Hertz  Pressure  and  Speed 
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Fig.  10  Low-Speed  Low-Temperature  MIL- L- 7808  Traction  Data 


-15- 


XI!-- 1400 


mmmmmm  Data 

mm  mm  mm  Theory 

Speed"8 !  448  in/8ec 
Temperature:  140‘F 

s'' 

• 

/ 

s 

/ 

/ 

/  ^ 

— - 

R 

• 

3,000  psi 

/ 

/ 

/  / 

// 
i  I 

s''  1 

PHZ  -  13 

0,000  psi 

tf 

// 

if 

y  / 

_ JU 

- ^ 

,/ 

s 

'^1 

_ 

- 

/ 

a 

s'' 
s  ^ 

PHZ-12 

3,000  psi 

/ 

^  * 

PHZ  *  10( 

3,000  psi 

j^s 

25  50  75  100  125  150 

Slip  Rate,  ips 


Fig.  11  Low-Speed  High-Temperature  MIL-L-7808  Traction  Data 
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Fig.  12  Medium-Speed  Low- Temperature  MIL-L-7808  Traction  Data 
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Fig.  13  Medium-Speed  High-Temperature  MIL- L- 7808  Traction  Data 
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Fig.  14  tilgh-Speed  Low- Temperature  MIL-L-7808  Traction  Data 
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The  base  viscosity  y0  obtained  from  MTI  viscosity  measurements  has  been 
fitted  with  a  viscosity  index  formula  of  the  form 

ViQ  *  1.45  x  10  ^{exp[exp(-3.705£n  T  +  460  +  24.39}]  -  .871  (3) 

2 

where  the  viscosity  y  is  in  units  lb-sec/in  .  This  formula  corresponds 
to  the  solid  line  passing  through  the  viscosity  data  appearing  in 
Figure  16. 

The  viscosity  data  shown  in  Figure  16  were  obtained  at  MTI  for  the 
MIL-L-7808  lubricant  used  in  all  of  the  traction  and  film  thickness 
studies.  This  lubricant,  which  was  supplied  by  the  Air  Force,  was 
identified  by  code  number  0-67-23.  Viscosity-pressure  data  was  ob¬ 
tained  by  Midwest  Research  Institute  (Reference  4)  for  a  different 
batch  of  the  same  lubricant  identified  by  Air  Force  code  number  0-67-20. 

The  MWRI  data  are  shown  in  Figure  17.  The  pressure  coefficients  of 
viscosity  were  obtained  from  the  slopes  of  the  £ny  versus  p  curves  in 
the  low-pressure  region  0  <  p  <  50,000  psi  as  shown  by  the  solid  lines 
in  Figure  17.  These  slopes  follow  the  relationship  with  temperature 

'  ,„-6  .  .013332  ,,x 

a  =  5.227  x  10  +  (4) 

where  a  has  units  in  /lb  and  T  is  in  °F. 

Since  film  thickness  is  largely  determined  in  the  inlet  zone  where  the 
pressures  are  somewhat  lower  than  in  the  center  of  contact  and  with 
pressure  coefficient  of  viscosity,  this  method  of  estimating  a  should 
be  adequate.  Equations  (2)  through  (4),  together  with  a  thermal  con¬ 
ductivity  of  .056  Btu/hr/f t/°F,  completely  describe  all  of  the  lubri¬ 
cant  property  input  used  in  obtaining  film  thickness  predictions. 

Traction  Predictions 

The  prediction  of  lubricant  traction  is  sensitive  to  lubricant  rheology 
at  the  high  pressures  occurring  at  the  center  of  contact.  The  low 
pressure  viscosity  data  used  in  determining  film  thickness  is  thus 
totally  inadequate  in  determining  traction.  Comparisons  between  trac¬ 
tion  for  5P4E  polyphenyl  ether  predicted  on  a  basis  of  low  pressure 
viscosity  data  with  measurements  were  performed  in  Reference  1  and 
showed  order  of  magnitude  discrepancies  between  theory  and  experiment. 
Similar  comparisons  have  been  made  here  for  MIL-L-7808  oil  with  simi¬ 
lar  results. 

In  order  to  correlate  the  traction  data,  it  was  necessary  to  extrapolate 
viscosity-pressure  data  outside  of  the  range  in  which  it  had  been  ob¬ 
tained.  The  hypothetical  viscosity  relationship  that  was  used  at  895  ips 
is  shown  graphically  by  the  broken  high-pressure  lines  in  Figure  17. 


These  high-pressure  lines  are  given  by  the  formula 


u 


T 


1,45  x  10  7  x(3)  exp  (-9,95  + 
+  ,36  x  10  ^  p) 


4,189  x  105 
(T  +  115)2 


(5) 


The  above  formula  represents  a  comprise  in  regard  to  reasonability  with 
respect  to  the  MWRI  viscosity  data  and  the  ability  to  obtain  a  reason¬ 
able  fit  to  the  traction  data. 

The  quantity  x(U)  is  a  rolling  speed  parameter  that  has  been  used  as  a 
multiplier  to  allow  for  the  observed  decrease  in  traction  with  increas¬ 
ing  rolling  speed  that  is  not  accountable  to  predicted  film  thickness 
variations  alone.  Equation  (5)  was  developed  by  taking  x(U)  to  be  1 
at  895  ips. 

Tractions  are  predicted  by  solving  the  equations  for  the  shear  stress 
tx  given  in  Reference  5  and  integrating  these  shear  stresses  over  the 
elliptical  contact  zone  with  the  use  of  the  assumption  that  the  pres¬ 
sure  is  Hertzian  as  was  performed  in  Reference  1,  The  equations  for 
predicting  traction  as  derived  in  the  above  references  are  summarized 
below.  The  Newtonian  shear  stress  relationship  is  given  by 


jju  _  \ 

3z  ~  bT 


(6) 


where  uT  is  the  viscosity  relationship  given  by  Equation  (5);  u  is  the 
local  lubricant  velocity,  and  z  is  a  coordinate  variable  which  varies 
across  the  film  -  h/2  <  z  <  h/2.  The  energy  equation  in  accordance 
with  Reference  5  assumes  the  form 


\  <u  ■  »> 


There  are  three  constraints  on  Equations  (6)  and  (7)  which  are 

T  =  T  ,  u  »  U,  at  z  =  h/2  (8) 

o  2 

and 

u  =  U  at  z  -  0  (9) 

Equation  (8)  may  be  regarded  as  providing  the  boundary  conditions  on 
the  second  order  system  given  by  Equations  (6)  and  (7)s  and  Equation 
(9)  may  be  regarded  as  the  additional  constraint  for  determi.— ;.g  t  „ 

Equations  (4)  through  (9)  together  with  the  Hertzian  pressure  assump¬ 
tion 
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and  the  elastohydrodynamic  performance  code  prediction  of  film  thick¬ 
ness  provide  all  of  the  necessary  equations  for  calculating  the  shear 
stress  ix  at  any  point  in  the  contact  zone.  Solutions  to  the  above 
system  are  readily  obtained  numerically  with  the  use  of  the  fourth 
order  Runge-Kutta  method. 

The  tractive  force,  F  ,  is  obtained  by  integrating  the  shear  stress 
over  the  Hertzian  contact  area  as  shown  mathematically  by 


ff 

n  J 


/l  -  (y/a)‘ 


t  dx  dy 
x 


which  may  be  readily  reduced  to  a  single  integral  by  expressing  the 
coordinates  in  polar  form. 

The  agreement  between  calculated  and  measured  tractions  at  the  extreme 
temperatures  at  895  ips  is  shown  in  Figures  12  and  13.  It  can  be  seen 
that  good  agreement  is  obtained  over  a  range  of  loads  and  temperatures 

The  quantity  x(U)  which  has  been  used  to  modify  viscosity  to  account 
for  possible  time  effects  is  shown  graphically  in  Figure  18.  The 
points  shown  are  those  actually  used  in  obtaining  the  theoretical 
curves  shown  in  Figures  10  through  15.  The  solid  line  represents  the 
formula 


X(U)  = 


(900)0» 23 

U 


U  900  ips 


U  >  900  ips 


Hence,  the  traction  correlation  may  be  extended  over  a  range  of  roll¬ 
ing  speed  as  well  as  temperatures  and  loads. 

The  above  correlation  is  not  intended  to  represent  a  full  theoretical 
description  of  lubricant  traction,  and  the  quasi-empirical  nature  of 
the  correlation  must  be  stressed. 

Complete  high-pressure-viscosity  data  for  the  fluid  tested  are  lack¬ 
ing  and  whereas  the  description  given  by  Equation  (6)  is  a  possible 
one,  it  does  not  necessarily  represent  true  viscosity  behavior. 
Similarly,  tractions  have  been  computed  with  predicted  film  thick¬ 
nesses  which  could  also  be  in  considerable  error.  In  light  of  these 


and  other  uncertainties,  the  analyses  used  here  should  be  regarded 
primarily  as  a  useful  correlating  device  for  interpolating  and  extra¬ 
polating  traction  data  over  reasonable  ranges;  but  care  should  be 
taken  before  drawing  any  theoretical  conclusions  regarding  the  high 
pressure,  high  shear  rate,  short-time  rheologica]  behavior  of  the  lub¬ 
ricant. 
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SECTION  III 

CAPACITANCE  FILM  THICKNESS  MEASUREMENT  AND  ANALYSIS 


Capacitive  film  measuring  techniques  have  been  described  in  the  literature  since 
1955-  Crook  (6),  during  the  late  fifties  and  early  sixties,  found  through  the  use 
of  the  capacitive  approach  many  of  the  unanticipated  film  thickness  relation¬ 
ships  between  rolling  speed  and  contact  loading. 

The  measurement  of  film  thicknesses  with  capacitive  techniques  requires  that  the 
test  roller  specimen  be  isolated  electrically  from  the  mating  roller.  The  me¬ 
thod  used  in  the  present  study  was  explained  in  detail  in  Reference  1.  The 
lower  disc  of  the  test  rig  was  designed  to  use  high  compression  plexiglass  as 
inserts  for  electrical  isolation.  The  choice  of  thickness  and  material  was 
governed  primarily  by  the  necessity  to  minimize  the  fixed  capacitance  to  ground. 
Direct  electrical  connection  to  each  roller  was  made  via  two  rotating  cables 
ending  in  two  mercury  baths  for  positive  electrical  integrity  to  the  capacitance 
bridge  system.  Any  attempts  to  use  electrical  paths  of  grounding  through  the 
shaft  support  bearings  resulted  in  an  erratic,  undesirable  ground  loop. 

The  total  capacitance  of  such  a  design  is  the  sum  of  the  isolated  disc  to  ground 
and  the  contact  zone  capacitance  formed  by  the  thin  oil  film.  Since  these  capa¬ 
citances  are  in  parallel 

C  =  c  +C 

total  constant  to  ground  contact  zone 

In  the  present  design  the  capacitance  of  the  contact  zone  was  equal  to  or 
greater  than  the  constant  capacitance  to  ground  when  the  test  roller  was  running 
under  loaded  conditions. 


The  exact  evaluation  of  the  capacitance  present  at  the  lubricated  contact  point 
of  the  discs  depends  on  a  knowledge  of  the  mating  disc  curvatures  and  the  way 
in  which  the  dielectric  coefficient  of  the  separating  material  varies  with  pres¬ 
sure  and  temperature.  Thus, 


C 


contact  zone 


C(K,h,T,PHZ) 


Several  independent  experiments  were  performed  to  check  out  the  various  indepen¬ 
dent  parameters  and  the  theoretical  expectations.  The  contact  zone  capacitance 
was  investigated  for  sensitivity  of  gap  separation  with  and  without  contact 
loading.  Lubricant  dielectric  coefficient  sensitivity  with  temperature  was 
also  investigaged . 

The  calculated  capacitance  of  the  contact  area  is  obtained  numerically  from 


C  „  ,  -  K  e  \  7  =  K  «  0 

contact  zone  o  J  h  o  h 

where  (A/h)  is  evaluated  numerically  with  the  ut,e  of  an  analysis  which  corrects 
for  Hertzian  elastic  distortion  both  interior  and  exterior  to  the  contact 
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ellipse  (see  Appendix  III  for  analyses).  The  integration  is  carried  to  in¬ 
finity  in  the  rolling  direction  but  truncated  at  the  disc  edge  perpendicular 
to  the  direction  of  rolling  (see  Figure  19).  Capacitive  contributions  from 
surfaces  beyond  the  edge  are  taken  as  a  constant  portion  of  the  fixed  isolated 
disc  to  ground  capacitance. 

1.  EXPERIMENTAL  SET-UP  AND  PRELIMINARY  CHECK  WITH  ANALYSIS 

Two  conditions  of  filling  the  gap  between  the  discs  with  dielectric  material 
were  considered  analytically.  The  first  case  considered  a  single  dielectric 
in  and  around  the  disc  whereas  the  second  considered  a  double  dielectric  fill¬ 
ing  of  the  gap  (see  Figure  20). 

A  preliminary  check  of  analytically-determined  capacitances  and  experimental 
variations  with  disc  separation  was  performed  with  the  test  arrangement  shown 
in  Figure  21.  The  contact  zone  gap  contained  a  single  dielectric  and  varied 
from  zero  to  one-hundred  mils  with  increments  as  small  as  2.8  microinch  in  the 
10  to  100  microinch  range  of  separation.  Gap  reading  accuracy  was  2.8  micro- 
inches  from  the  standard  dial  micrometer  bed  and  lever  setup  shown.  The  lower 
shaft  was  isolated  electrically  from  the  grounded  upper  shaft. 

The  contact  zone  capacitance  as  determined  experimentally  and  analytically  as 
a  function  of  minimum  gap  separation  is  plotted  in  Figure  22.  The  fixed  capa¬ 
citance  to  ground  was  taken  as  the  capacitance  of  the  system  determined  with  a 
gap  of  100  mils  and  was  subtracted  from  each  subsequent  reading.  The  measured 
capacitances  were  obtained  with  a  standard  carrier  frequency  capacitance  bridge. 
The  auto  balance  precision  bridge  combines  speed  of  measurement  with  .01  per¬ 
cent  accuracy  and  a  digital  display  of  capacitance  to  five  significant  places, 
as  well  as  the  conductance  of  the  component  under  test.  The  bridge  used  has 
seven  decade  ranges  from  micro- farads  to  pico- farads.  Each  decade  provided 
digital  capacitance  readings  or  analog  signals  with  slightly  less  accuracy  for 
automatic  graphic  display  of  measured  capacitances. 

In  addition  to  knowing  the  unloaded  capacitive  dependence  on  separation,  the 
fully  loaded  sensitivity  must  also  be  determined  before  lubricant  film  thick¬ 
ness  can  be  found  from  capacitance  data.  Contact  zone  capacitance  is  very  sen¬ 
sitive  to  loading  since  the  area  of  the  Hertzian  zone  increases  as 

2/3 

Contact  Area  ot  Load 

An  experimental  check  of  the  expected  capacitive  changes  with  increasing  loads 
was  performed  with  mica  separating  shims.  Ihin  cleavages  1.3  mils  thick  were 
successively  stacked  between  the  two  stationary  discs.  Increases  in  capacit¬ 
ance  with  higher  loads  and  decreases  in  capacitance  with  a  greater  number  of 
mica  shims  were  observed.  Tne  results  are  displayed  in  Figure  23.  The  use  of 
mica  as  a  load  calibrating  material  presents  some  problems  experimentally.  The 
dielectric  coefficient  can  vary  depending  upon  the  trace  mineral  content  of 
the  pieces  used  in  testing.  Published  variations  from  5.0  to  8.0  are  common 
and  are  not  accurate  enough  for  calculating  purposes.  The  two  points  at  80  ksi 
indicated  in  Figure  23  were  used  in  determining  ground  capacitance  results 
and  the  dielectric  coefficient  of  the  mica.  The  dielectric  coefficient  was  de¬ 
termined  to  be  7.45  which  lies  within  the  range  specified  above.  The  coales¬ 
cence  of  theory  and  experiment  for  both  loaded  and  unloaded  separations  of  the 


Uniformly  Thick  Rigid 
Non-Conductive  Shim 


Fig.  20  Dielectric  Conditions  Considered  Analytically 
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test  specimens  provided  enough  evidence  to  proceed  in  measuring  full  fluid 
film  thicknesses  under  rolling  conditions. 

2.  ANALYTICAL  ASSUMPTIONS 

Three  basic  assumptions  were  made  in  making  the  film  thickness  determinations 
from  the  capacitance  data.  Although  full  Hertzian  distortions  were  used  in 
calculating  the  contact  area  and  gap  distribution,  the  trailing  edge  dip  pre¬ 
sent  under  rolling  conditions  was  neglected.  Since  the  dip  was  observed  in 
the  optical  work  (see  Section  IV)  to  extend  over  a  minor  portion  of  the  con¬ 
tact  area,  this  would  result  in  a  small  percentage  inaccuracy  in  the  relative 
capacitances  obtained.  A  twenty  percent  trailing  edge  dip  occupying  ten  per¬ 
cent  of  the  contact  zone  area  could  be  expected  to  have  no  more  than  two  per¬ 
cent  error  in  expected  capacitance  produced  by  the  contact  zone. 

The  second  simplifying  assumption  is  related  to  the  distribution  of  the  loci 
of  all  points  of  equal  separation  in  and  around  the  contact  zone  between  the 
discs  referred  to  here  as  "fringp  lines".  Feints  of  equal  separation  are  cal¬ 
culated  along  the  major  and  minor  axes  in  and  around  the  load- carrying  zone  of 
the  test  specimens i  Although  the  fringe  lines  are  not  truly  elliptical 
(a  maximum  radial  difference  of  approximately  6  percent  was  observed  analyti¬ 
cally),  they  were  assumed  to  be.  For  example  ,  the  numerical  point  by  point 
distributions  are  compared  to  their  elliptical  counterparts  and  shown  in  the 
schematic  of  Figure  24.  Elliptical  distributions  were  assumed  only  to  reduce 
calculating  machine  time.  No  significant  error  in  calculating  is  generated 
with  a  maximum  disparity  of  6  percent,  since  the  total  capacitance  is  a  summa¬ 
tion  of  all  differential  areas  with  very  small  curvatures. 

The  third  conceptual  simplification  made  in  calculating  capacitances  is  related 
to  the  temperature  within  the  rolling  contact.  The  rolling  zone  was  considered 
to  have  a  sufficiently  small  temperature  rise  such  thar  it  could  be  neglected 
for  the  loads  investigated.  Ibis  assumption  is  valid  since  the  temperature 
variations  in  the  dielectric  coefficient  were  found  to  be  smaller  than  .04  per¬ 
cent  per  F  of  a  temperature  rise  from  70°F  t-°  300°F  for  the  MlL-L-  7808  tested. 
It  has  also  been  observed  by  Hamilton  and  Moore  (7  ;  that  tr.e  temperature  rise 
in  a  pure  rolling  contact  similar  to  tne  present  setup  is  cf  the  order  of  1°F. 
Greater  temperature  rises  would  occur  in  sliding  contacts,  however,  even  100°F 
increases  in  temperature  would  not  significantly  affect  tne  results.  In  fact, 
experimentally  some  tests  were  observed  with  and  without  slip  up  to  10  percert 
and  very  little  change  in  capacitance  (~  2  percent)  was  observed. 

3.  DIELECTRIC  PRESSURIZATION  AND  CALIBRATION  CURVES 

Dielectric  coefficients  of  non-polar  fluids  are  known  to  respond  to  mass  density 
changes  as  a  result  of  applied  pressure.  The  Clausius-Mossotti  (8)  equation  is 
normally  assumed  in  describing  the  reactive  changes  unless  precise  pressure- 
dielectric  measurements  are  made  on  the  test  lubricant.  Direct  measurements 
by  Galvin,  Naylor  and  Wilson  (9)  have  produced  an  explicit  relation 

K/K  =  1  +  2.59  x  10  6  p  •  29.8  x  10"12  p2 

O  r  r 

for  a  mineral  oil  where  p  is  the  pressure  applied  to  the  fluid  in  pounds  per 
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square  inch.  This  estimation  of  the  effective  point  wise  dielectric  coeffi¬ 
cient  throughout  the  contact  region  was  assumed  in  the  present  film  thickness 
determinations.  It  is  apparent  that  pressure  levels  up  to  80,000  psi  would 
not  significantly  effect  the  dielectric  coefficient.  Film  thicknesses  below 
100  microinches,  as  determined  with  and  without  the  above  pressure  modifica¬ 
tions,  are  presented  in  Figure  25  for  comparison.  Shown  are  two  analytically- 
derived  sets  of  calibration  data  from  which  film  thicknesses  can  be  taken  for 
any  level  of  disc  loading.  The  dotted  curves  contain  all  significant  load 
distortions  and  separation  effects  but  do  not  include  the  point  by  point  pres¬ 
surization  effect  on  the  lubricant  dielectric  coefficient.  The  solid  set  of 
curves  include  the  dielectric  pressurization  effect  for  each  load  depicted. 

..^though  it  appeared  at  first  there  was  a  possibility  of  getting  a  signficant 
change  when  the  dielectric  pressurization  was  considered,  the  calculations  are 
less  sensitive  than  was  anticipated  for  loading  up  to  the  170,000  psi  range. 

Since  there  was  shown  with  air  to  be  a  valid  agreement  between  experimentally 
determined  and  calculated  films  for  disc  separations  down  to  30  microinches, 
five  points  between  1  and  10  mils  disc  separation  were  used  for  defining  the 
no  load  under  fully  flooded  conditions  on  the  disc  machine.  A  full  set  of 
analytically  obtained  capacitance  data  and  the  corresponding  range  of  thick¬ 
nesses  for  which  experimental  values  were  obtained  is  presented  in  Figure  26. 

A  similar  set  of  curves  obtained  analytically  for  polyphenyl  ether  is  dis¬ 
played  in  Figure  27. 

4.  TEST  DATA 

The  capacitance  data  as  taken  from  the  disc  machine  is  plotted  automatically 
from  analog  signals  as  a  function  of  rolling  speed.  A  typical  plot  of  actus* 
recorded  capacitance  data  is  displayed  in  Figure  28.  For  the  load  and  speed 
indicated,  the  MIL-L«7808  lubricant  plot  is  not  valid  at  disc  rotational  fre¬ 
quencies  less  than  800  rpm.  The  film  at  this  point  has  become  sufficiently 
small  to  allow  some  degree  of  asperity  coitacting  which  prevents  the  capacit¬ 
ance  bridge  from  operating  properly.  Refer  to  Section  VI  for  plots  which  show 
the  nature  of  the  asperity  contact  interactions  below  1000  rpm.  The  capacit¬ 
ance  obtained  at  each  rolling  speed  is  then  used  analytically  to  get  true  film 
thickness  as  a  function  of  disc  load  and  speed.  Successive  runs  made  at  ele¬ 
vated  temperatures  show  the  dependence  of  film  thickness  on  fluid  viscosity. 

Loading  effects  on  MIL-L-7808  film  thickness  have  been  observed  by  various 
experimenters  in  the  past.  Kannel,  Bell  and  Allen  (3  )  presented  evidence  for 
a  large  load  effect.  The  loading  sensitivity  in  the  present  study  made  with 
capacitive  techniques  is  functionally  dependent  upon  the  load  level  under  ex- 
amir  ition.  At  heavy  loads  the  film  thickness  was  found  to  vary  with  load  to 
the  .2  power  or  less  whereas  at  somewhat  lighter  loads  power  variations  as 
high  as  .6  were  encountered. 

From  a  series  of  capacitance  curves  (similar  to  the  one  shown  previously  in 
Figure  28)  taken  at  different  loads,  a  functional  plot  of  observed  film  thick¬ 
nesses  versus  rolling  speed  may  be  obtained  with  the  analysis.  The  observed 
film  versus  speed  relation  is  the  contents  of  Figure  29. 
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Fig.  26  Analytical  and  Calibration  Data  for  Capacitance  Versus  Film 
Thickness 
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Examination  of  the  slopes  of  the  various  film  thickness  vs.  speed  lines  ob¬ 
tained  at  78°F  from  the  MIL-L-7808  indicates  that  film  thickness  varies  with 
speed  to  the  0.79  power  at  80  ksi  and  decreases  to  a  0.65  power  variation  at 
140  ksi.  These  results  are  quite  comparable  with  the  0.7  put  forth  by  Dowson 
and  Higginscn  (10),  A  wider  i  s.i&e  of  power  variations  were  observed  from  the 
polyphenyl  e^hor  data.  The  observed  film  thickness  varied  with  speed  to  the 
.85  power  with  light  loads  and  .40  power  with  heavy  loads. 

Similar  recordings  of  capacitance  data  were  taken  with  polyphenyl  ether  as  a 
contact  lubricant.  A  comparison  of  film  thicknesses  obtained  with  polyphenyl 
ether  as  the  lubricating  media  show  the  test  MIL-L-7808  maintains  a  thinner 
film  for  every  load  and  speed  measured.  See  Figures  30  and  31.  In  fact,  at 
room  temperature  the  observed  film  thickness  for  polyphenyl  ether  is  nearly 
constant.  Capacitive  outputs  are  irratic  when  the  polyphenyl  was  above  200°F 
and  provide  little  film  thickness  information  above  that  temperature. 


Fig.  30  Cap 
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SECTION  IV 


OPTICAL  FILM  THICKNESS  MEASUREMENT  AND  ANALYSIS 


The  design  details  of  the  MTI  optical  test  rig  are  described  in  a  previous  year's 
annual  report  (1).  The  first  significant  lubricated  contact  experiments  made 
with  optics  were  performed  by  Cameron  (11).  New  approaches  (12)  in  recent  years 
have  upgraded  the  quality  of  optical  film  thickness  measurements  in  loaded  roll¬ 
ing  contacts.  The  advantage  of  the  photographic  technique  is  that  a  single  pic¬ 
ture  reveals  the  full  relative  topological  distribution  in  and  around  the  observed 
contact  region.  If  the  finishes  and  coatings  are  carefully  prepared,  measured 
variations  of  film  thickness  as  small  as  one-eighth  of  a  wave  length  of  light 
are  common  (typically  2  or  3  microinches). 

1.  TEST  SET-UP 

The  complete  test  assembly  was  an  intricate  arrangement  which  provided: 

(a)  Visual  observance  of  the  elliptical  contact  region 

(b)  Monochromatic  filtered  illumination  of  the  contact  region 

(c)  Variable  speed  control  and  detection  of  the  rolling  disc  specimens 

(d)  Variable  loading  capability  of  mating  test  discs 

(e)  Stroboscopic  illumination  for  photographic  recording  of  the 
interference  fringe  pattern 

(f)  An  electrical  triggering  device  to  synchronize  the  optical  access 
port  with  the  illuminating  strobe  light. 

The  MTI  optical  setup  used  curved  surface  geometries  for  both  the  metal  and  trans¬ 
parent  quartz  mating  discs.  Before  the  present  study  only  limited  optical  data 
related  to  lubrication  experimentation  have  been  available  for  geometries  other 
than  steel  balls  on  an  optical  flat  arrangement.  The  optical  quality  trans¬ 
parent  disc  was  an  uncrowned  cylinder  with  a  three  inch  outside  diameter  and  a 
two  inch  inside  diameter  (see  Reference  1  for  photographs  of  specimens).  The 
quartz  disc  was  mounted  with  epoxy  onto  the  disc  machine  lower  shaft  assembly 
as  shown  in  the  cross-sectional  diagram  of  Figure  32-  The  shaft  to  which  the 
transparent  cylinder  was  attached  was  hollow  for  optical  access  at  three  points 
around  the  disc  periphery. 

The  interference  technique  used  in  the  present  study  makes  use  of  thin  film 
coating  techniques  with  chromium  in  order  to  increase  contrast  of  the  interfer¬ 
ence  fringes.  A  20  percent  reflective  coat  was  deposited  upon  the  outside 
diameter  of  the  quartz  disc.  Standard  microscopic  components  which  were  opti¬ 
cally  coated  for  reducing  excess  internal  reflections  were  used  in  the  main 
body  of  the  optical  setup. 
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Transparent  Quartz  Disc 


Loading  was  accomplished  through  the  shaft  bearings  which  supported  each  test 
specimen.  This  eliminated  the  common  lubricant  heating  problem  which  arises 
when  a  sliding  ball  and  cup  made  of  teflon  is  used  for  contact  loading  purposes. 

The  test  lubricants  examined  were  the  same  as  used  in  obtaining  the  traction  and 
capacitive  data  described  in  the  previous  sections.  The  physical  properties 
related  to  the  test  lubricant  are  given  for  reference  in  Appendix  I.  The  syn¬ 
thetic  turbine  oil  of  MIL-L-7808  specification  provided  the  best  optical  data 
of  the  two  test  fluids  and  allowed  a  comparison  of  nominal  film  thickness 
obtained  by  the  separate  methods  of  capacitive  and  optical  measurement.  Attempts 
to  produce  fringes  with  polyphenyl  ether  yielded  unsatisfactory  results.  This 
was  most  likely  related  to  the  speed  range  with  which  the  test  rig  could  be 
varied.  The  polyphenyl  at  room  temperature  has  a  viscosity  of  approximately  600 
centistokes  and  provided  extremely  high  film  thickness  {>  100  microinch)  at  the 
lowest  rolling  speed  obtainable  on  the  present  test  rig. 

A  schematic  of  the  illumination  technique  used  in  the  test. rig. is  displayed  in 
Figure  33.  The  strobe  unit  used  for  photographic  images  was  a  single  pulse 
unit  with  1/10,000  second  duration  time-  The  unit  was  triggered  via  the  "bulb" 
shutter  position  on  the  camera  and  a  pair  of  sliding  pin  contacts  on  the  shaft 
which  held  the  quartz  specimen.  A  yellow  monochromic  transmission  filter  was 
used  between  the  light  source  and  the  contact  zone.  The  transmitted  yellow 
light  wave  length  was  independently  measured  with  a  dual  slit  diffraction  test 
and  the  slit  spacing  determined  with  the  6328A  light  of  a  helium-neon  laser. 

The  predominate  wave  length  of  the  contact  illumination  source  after  filtering 
was  found  to  be  5890A  <23.2  microinches)  in  air. 

Visual  observations  of  the  contact  zone  were  also  made  with  a  rapid  cycle  multi¬ 
flash  unit  of  lower  light  intensity  output.  That  system  was  also  triggered 
synchronously  with  the  shaft  and  light  port  alignment  allowed  visual  observation 
of  the  contact  interference  pattern  as  it  changed  while  adjusting  the  rolling 
speed. 

The  drive  system  was  synchronized  during  the  optical  work  for  pure  rolling  op¬ 
eration.  No  investigations  with  slip  between  the  discs  were  made  Shaft  speeds 
were  monitored  with  a  Fotonic  light  probe,  a  strjped  decoding  band,  and  a  four- 
place  digital  counter.  Rotational  frequencies  as  low  as  10  rpm  could  be  maintained. 
However,  motor  cogging  of  the  variable  speed  unit  at  such  slow  rotational  fre¬ 
quencies  limited  the  useful  speed  range  to  20  rpm  or  greater. 

2.  OPTICAL  DISTORTION  CORRECTION  OF  LLLIF7IC4L  BANDS 

The  index  of  refraction  used  throughout  the  experiment  were: 

Material  Refractive  Index 

Polyphenyl  Ether  (SP'-.F,  1.630 

MIL-L-7808  1-395 

Quartz  Test  Disc  1.458 

and  was  assumed  constant.  The  quartz  disc  has  a  modulus  of  elasticity  of 
approxfmately  one-half  of  that  of  the  steel  disc  and  therefore  requires  Hertz 
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stresses  of  less  than  50,000  psi  for  equivalent  contact  distortions  equal  to 
the  100,000  psi  level  with  two  steel  mating  discs.  In  fact,  the  choices  of 
loading  were  made  such  that  the  minor  axis  (along  the  rolling  direction)  of  the 
contact  ellipse  were  the  same  for  some  tests  made  with  the  capacitive  technique. 


Hertz  stresses  of  50,000  psi  would  result  in  less  than  a  4  percent  change  in 
the  refractive  index  and  was  not  considered  significant- to  adjust  any  film 
thickness  data  for  pressurization.  The  4  percent  change- would  only  occur  over 
a  limited  area  of  the  contact  zone  and  was  calculated  from  the  classic  Lorenz- 
Lorenz  equation 

,1  +  2A.1/2 
n  =  ( - _) 

P  1  +  A 

where, 


n  2  -  1 


A  =  —  (— « 


o  n  +2 
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with  measured  densities  taken  from  the  Midwest  data  in  Appendix  I. 
sures  up  to  50,000  indicate 


Hertz  pres- 


n  -  n 

J - °  <  4% 

n  — 

o 

Because  of  the  geometry  of  the  transparent  test  disc,  it  became  necessary  to 
apply  a  constant  distortion  correction  factor  to  each  observed  fringe  position 
along  the  contact  rolling  direction  in  order  to  interpret  the  actual  spacing. 
The  spacing  perpendicular  to  the  rolling  direction  needed  no  correction.  The 
image  distortion  caused  by  using,  a  cylindrical  roller  instead  of  the  optical 
flat  used  by  previous  investigators  is  derived  in  Appendix  IV.  The  correction 
factor 

f  _  Observed  Fringe  Position 
D  “  Actual  Fringe  Position 

is  found  from  the  optical  geometry  to  be 

nG 

f°  ^.(nG  -  1)  +  1 
r  G 

The  disparity  between  observed  and  actual  dimensions  is  schematically  represen¬ 
ted  in  Figure  34.  An  understanding  of  the  distortion  expected  can  be  obtained 
from  the  side  view  of  the  test  discs  shown  in  Figure  35.  The  distortion  factor 
was  evaluated  to  be 


fD  =  .8641 

for  the  optical  geometry  or  approximately  14  percent  between  observed  and  actual 
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spacings.  This  value  was  used  to  correct  all  observed  dimensions  in  the  rolling 
direction  and  is  inherent  in  all  photographs  shown  in  this  report. 

In  order  to  test  the  optical  equipment  and  analysis,  a  preliminary  static  load¬ 
ing  test  was  performed  with  MIL-L-7808  in  the  contact  zone  of  the  mating  discs. 
Six  loads  giving  Hertz  pressures  from  22,000  psi  to  49,000  psi  provided  photos 
similar  to  those  shown  in  Figure  36.  The  increase  in  contact  zone  area  is 
apparent  with  increased  loading  although  the  actual  points  of  metal  to  quartz 
contact  within  that  zone  cannot  be  seen  in  the  photographs.  Exhibited  are  the 
fringes  of  equal  gap  separation  between  the  loaded  discs  for  gaps  greater  than 
one-eighth  of  a  wave  length  of  the  illuminating  yellow  light. 

The  interference  intensity  at  the  chromium  coat  depends  upon  the 

(a)  Chromium  film  thickness  deposited 

(b)  Relative  index  of  refraction  compared  to  that  of  quartz 

(c)  Separation  gap  between  the  discs 

(d)  Phase  change  of  the  incident  light  at  the  polished  metal  surface. 

The  quartz  mating  disc  used  was  coated  with  a  20  oercent  reflective  coating  of 
chromium  which  increases  the  contrast  of  the  light  interference  bands  in  the 
contact  region.  The  chosen  coating  was  assumed  to  induce  an  even  number  of 
it  phase  shifts  in  the  reflected  light  as  indicated  by  Holden  (13)  from  his  thin 
film  work  with  silver  coatings.  This  means  that  the  first  surface  reflection 
has  the  same  phase  as  the  incoming  light  since  the  index  of  refraction  of  the 
oil  is  less  than  that  of  the  quartz  used.  The  reflected  light  from  the  polished 
crown  metal  disc,  however,  undergoes  a  phase  shift  equal  to  (  it  )  according  to 
Andrews  (14)  which  is  the  equivalent  of  one-half  of  a  wave  length  of  light. 

Light  which  has  traveled  a  gap,  h,  such  that 


where, 

i  -  0,  1,  2,  3 

wlienl  is  t! Ewave  length  of  light  in  the  separating  lubricant,  causes  bright  fringes 
at  the  chromium  coat  whenever  the  index  of  refraction  of  the  lubricant  is  less 
than  the  refractive  index  of  the  transparent  disc.  In  the  present  work  dark 
fringes  appear  whenever  the  gap,  h,  is  such  that 

h  -  (|i)  J 

where, 

i  -  1,  2,  3 . 

and  X  is  again  the  wave  length  of  the  illuminating  light  in  the  test  lubricant. 


-54- 


The  inner  most  dark  band  which  occurs  is  a  special  case  of  interference.  There 
is  enough  lubricant  trapped  between  the  discs  to  allow  a  half  wave  phase  shift 
in  the  light  reflected  at  the  chromium  coat  but  not  enough  separation  to  provide 
a  half  wave  length  travel.  Therefore,  an  interference  minimum  or  a  dark  band 
appears.  The  central  contact  region  contains  little  or  no  lubricant.  The 
metal  is  so  close  to  the  chromium  coat  that  very  little  light  is  transmitted 
and  all  light  reflected  is  in  phase  so  a  central  maximum  or  light  gray  area 
appears.  A  summary  of  the  variations  of  intensity  expected  from  a  thin  film 
gap  is  given  by  Figure  37. 

The  inner  moot  bright  ring  (shown  previously  in  Figure  36)  occurs  when  the 
separation,  h,  between  the  discs  is 


h  = 


air 


4n 


oil 


=4.15  microinches. 


3.  TEST  DATA  WITH  ANALYTICAL  EXPECTATIONS 

The  actual  lo^d  carrying  contact  area  (such  that  h  =  0)  where  micro-asperities 
support  the  normal  load  between  the  discs  is  well  within  the  inner  most  white 
band.  The  analytically-calculated  dimensions  of  the  major  and  minor  axes  of 
the  contact  zone  separation  of  4.15  raicroinch  points  for  six  lead  levels  are 
compared  in  Figure  38.  The  results  of  that  figure  indicate  that  the  elasticity 
analysis  used  in  describing  the  distribution  of  surface  gaps  in  and  around  the 
contact  area  (given  in  Appendix  III)  are  sufficient  to  accurately  define  the 
experimentally  observed  distributions  with  load  under  static  contact.  In  fact, 
the  locus  of  all  points  of  equal  separation  for  several  loads  were  numerically 
evaluated  and  computer  plotted.  The  results  along  with  their  corresponding 
contact  photographs  are  shown  in  Figures  39  and  40.  Overlaying  the  computer 
plots  on  the  corresponding  photographs  reveal  a  positive  matchup  between  theory 
and  experiment. 

Over  three  hundred  photographs  were  taken  of  the  contact  region  under  various 
load-  speed  conditions.  The  photographs  shown  in  Figures  41  and  42  reveal  how 
the  film  pattern  changes  with  increasing  rolling  speed.  As  the  speed  of  the 
disc  increases,  the  nominal  film  increases.  The  central  minimum  film  region 
responds  by  changing  alternately  from  light  to  dark  and  back  again  as  the  speed 
is  brought  to  higher  and  higher  levels.  As  shown  in  Figure  43,  the  outer  bands 
move  toward  the  center  as  the  rolling  speed  increases.  What  was  a  dark  ring 
about  a  central  bright  region,  becomes  a  dark  central  contact  region  if  the 
nominal  gap  increases  by  a  quarter  of  a  wave  length  of  the  illuminating  light. 

A  photograph  showing  a  rolling  contact  taken  at  a  speed  where  a  dark  band  has 
moved  in  to  lie  partly  in  the  central  region  and  partly  along  the  trailing  edge 
dip  is  displayed  in  Figure  45c 

Photographs  reproduced  in  Figures  44  through  46  were  taken  as  the  disc  rolling 
Stoed  was  increased  from  26  rpm  to  488  rpm  with  a  Hertz  pressure  of  22,000  psi. 
A  •  000  psi  steel  on  quartz  pressure  was  equivalent  to  44,000  psi  for  the 

on-steel  test  arrangement. 

A  summary  of  a  second  set  of  variable  speed  tesrs  performed  with  a 

36,  .  Hertz  stress  is  presented  in  Figure  47.  The  choree  of  load  was  made 
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so  as  to  -<atch  the  contact  zone  minor  axis  dimension  in  the  rolling  direction 
to  the  oceel-on-steel  contact  dimensions  of  the  capacitive  test  run  at  80,000 
psi.  From  both  visual  and  photographic  observations  of  the  contact  region  over 
the  full  range  of  speeds  from  0  to  1,400  rpm  the  plot  shown  in  Figure  48  was 
compiled.  The  plot  depicts  five  variable  speed  tests  made  with  the  MIL-L-7808; 
three  capacitive  runs  and  two  optical  runs  under  separate  load  conditions. 

Although  the  film  thickness  speed  data  for  the  capacitive  observations  show  a 
marked  load  dependence  for  the  lighter  loads ,  this  is  not  evident  from  the 
optical  data.  The  optical  work  indicates  a  simple  film  speed  exponent  relation¬ 
ship  of  0.85  over  the  rolling  speed  range  of  30  rpm  (4.71  in/sec)  to  1,400  rpm 
(220  in/sec)  tested. 

The  analytical  line  shown  in  Figure  48  is  the  EHD  performance  code  prediction 
corresponding  to  the  conditions  under  which  the  36  kpsl  optical  measurements 
were  taken.  The  performance  code  prediction  reduces  very  nearly  to  the  Dowson 
and  Higginson  formula  under  these  conditions  since  inlet  zone  heating  effects 
are  predicted  to  be  small. 

The  optical  data  seems  to  substantiate  conventional  EHD  predictions  regarding 
load  sensitivity  at  low  pressures  whereas  the  capacitance  data  obtained  at 
higher  pressures  and  higher  speeds  seems  to  indicate  a  more  significant  load 
dependence. 
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SECTION  V 

X-RAY  FILM  THICKNESS  MEASUREMENT 


R 


The  X-ray  film  thickness  measuring  technique  was  initially  developed  by  workers 
at  Battelle  Memorial  Institute  in  the  early  1960 's  The  work  has  since  been 
heralded  as  the  first  valid  proof  of  the  existence  of  very  thin  films  in  roll¬ 
ing  element  contacts.  The  rolling  disc  machine  described  in  a  previous  Air 
Force  report  (Reference  1)  was  designed  to  allow  similar  X-ray  measurements 
The  conceptual  layout  is  diagrammed  in  Figure  49  X-rays  from  a  natural 
decay  source  are  passed  through  the  lubricating  film  of  the  disc  machine  in 
the  direction  of  rolling 

A  beam  of  X-rays  with  a  width  of  .090  inch  at  the  contact  zoi,~  was  passed  be¬ 
tween  the  .-oilers  and  detected  with  a  photomul » ip  1 y ing  tube.  The  dimensional 
characteristics  and  beam  arrangement  are  shown  in  Figures  50  and  51. 

The  radioactive  source  was  chosen  for  its  desirable  X-ray  energy  characteris¬ 
tics,  Figure  52  displays  the  energy  spectrum  emitted  from  the  Cadmium  109 
decay  source-  The  predominate  X-ray  peak  appears  at  22  kev  and  is  emitted  from 
the  metastable  silver  state  produced  on  decay.  The  22  kev  X-ray  wave  length  is 
very  close  to  that  of  the  molybdenum  target  X-rays  generated  for  the  pioneering 
work  originally  performed  by  Sibley,  Bell,  Orcutt,  and  Allen  (151  The  radia¬ 
tion  of  the  molybdenum  target  used  by  them  exhibits  peaks  fr am  the  Ka  and  Kg 
shells  which  have  energies  of  16  5  kev  and  18.6  kev  respectively-  X-ray  ener¬ 
gies  less  than  30  kev  are  desirable  since  they  penetrate  *-he  lubricants  quite 
well  and  yet  fill  not  penetrate  bearing  steels  to  any  appreciable  extent 

The  X-ray  detection  instrumentation  was  identical  to  the  original  sc  cup  used 
in  the  Battelle  work  A  schematic  of  the  electrical  linkup  used  for  monitoring 
and  calibrating  is  shown  in  Figure  53  ana  a  complete  list  of  the  test  equipment 
is  provided  in  Appendix  V  The  calibration  procedure  involved  the  determina¬ 
tion  of  absorption  coefficients  of  the  X-rays  in  the  test  lubricants,  as  well 
as  for  count  rate  determinations  for  fixed  disc  separations 

If  I  is  the  intensity  of  X-rays  on  a  piece  of  material  cf  thickness  dx,  and  dl 
is  the  decrease  in  intensity  of  the  X-rays  due  to  absorption  within  the  mater¬ 
ial,  then  u  the  absorption  coefficient  is  defined  such  that 

~  d  I 

4  “  Idx 

Thus , 

I  =  1  exp  (-  ux  , 


ib  che  .\-ray  intensity  vh-  a  >*  =  0  The  experimentally  determined 
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values  of  p.  found  from  the  two  test  fluids  and  X-ray  source  used  were 


Test  Material 


MIL- L- 7808 


Cd 


109 


X-Ray  Absorption  Coefficient 


-  .63  in 


Polyphenyl  Ither  (5F4E) 


.82  in 


The  calibration  of  the  disc  separation  was  performed  in  the  same  manner  as  with 
the  capacitance  techniques.  The  count  rate  of  X-rays  was  determined  for  fixed 
separations  between  the  discs  at  1,  2.  3,  5,  and  10  mils  separation.  Since  the 
beam  width  of  the  X-ray  source  was  confined  to  less  than  90  mils  both  in  height 
and  width  at  the  contact  zone,  the  calibration  was  made  with  no  load  on  the 
discs.  The  relationship  between  separation  and  count  rate  was  found  to  be  both 
linear  and  a  first  power  dependence  as  shown  by  the  log-log  plot  in  Figure  54. 

An  extrapolation  of  the  dependence  into  the  microinc'n  region  was  taken  for  the 
actual  determination  of  film  thickness  from  the  count  rate  data  taken  under  full 
load  and  rolling  speed  conditions. 

Well  over  150  sets  of  count  rate  data  under  dynamic  conditions  were  recorded 
during  testing.  At  least  five  recordings  of  X-ray  count  rate  were  taken  for 
each  of: 

(a)  Three  rolling  conditions:  zero,  470,  and  920  in/sec 

(b)  Five  loads:  100,  120.  140,  160,  and  172  kpsi 

(c)  Two  radiation  conditions:  beam-on  and  background  (beam-off) 

after  beam  alignment  had  been  achieved.  The  disc  machine  contact  zone  was 
located  with  metal  bluing  and  aligned  with  the  exit  hole  of  the  radiation  con¬ 
tainment  pig.  A  through  hole  m  the  container  was  opened  whenever  the  Cadmium 
109  source  was  rotated  to  the  OFF  position  in  the  container.  This  through 
hole  provided  a  convenient  way  (by  insertion  of  a  metal  rod)  to  physically 
align  the  X-ray  beam  both  in  the  vertical  and  horizontal  directions  with  the 
contact  zone. 

The  calibration  of  count  rate  with  disc  separation  was  determined  when  air 
filled  the  gap  between  the  discs.  The  actual  situation  under  rolling  is  similar 
to  that  shown  schematically  in  Figure  55,  where  the  X-rays  must  pass  through  a 
finite  thickness  of  the  test  lubricant.  The  precise  thickness  is  much  greater 
than  the  minor  axis  of  the  contact  zone,  however,  and  was  estimated  to  be 
s  .125  inches  from  visual  observation  for  the  flooded  conditions  of  actual 
running. 

Net  count  rate  data  used  to  determine  lubricant  film  thicknesses  were  obtained 
by  subtracting  the  background  recordings  taken  under  identical  loading  condi¬ 
tions  with  no  rocation  of  the  discs.  The  mean  values  of  film  thicknesses 
obtained  for  each  of  the  five  load  levels  investigated  are  plotted  in  Figure 
54  along  the  previously  discussed  calibration  line.  The  nominal  film  thick¬ 
nesses  for  rolling  speeds  of  920  in/sec  and  470  in/sec  were  found  to  be  ap¬ 
proximately  30  microinches  and  were  determined  from  the  two  mean  values  of  all 
load  data  of  the  separate  speeds  investigated.  Although  the  attenuation  of 
the  count,  rate  signal  was  know,  it  was  not  used  to  modify  the  data.  The  X-ray 


attenuation  is  exponential  with  traversed  length  in  lubricant  and  a  lubricant 
path  less  than  1/8  inch  would  not  significantly  change  the  observed  film 
thicknesses.  Any  consideration  of  the  loss  in  X-rays  in  the  lubricant,  how¬ 
ever,  would  be  interpreted  to  give  larger  films  than  the  30  microinches  found. 

No  correlation  between  variations  in  film  thickness  with  speed  and  load  could 
be  obtained.  Some  of  the  lower  loads  revealed  smaller  films  than  some  of  the 
higher  loads  investigated  and  the  mean  film  thickness  of  the  slow  speed  was 
larger  than  that  obtained  at  the  higher  test  speed.  The  partial  inconsisten¬ 
cy  in  the  data  was  assumed  to  be  a  result  of  a  low  signal-to-noise  ratio  ob¬ 
tained  electronically  from  the  available  X-ray  count  rates.  The  root-mean- 
square  of  the  data  gathered  under  all  loading  conditions  is  shown  also  in 
Figure  54  as  a  measure  of  the  scatter  observed  in  the  X-ray  data. 
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SECTION  VI 


ASPERITY  INTERACTION  AND 
PARTIAL  ELASTOHYDRODYNAMICS 


Observation  of  asperity  interaction  was  performed  wi th  the  aid  of  the  electronic 
"Asperitac"  developed  specifically  for  the  work  reported  here.  An  overview  of 
the  operational  characteristics  of  the  electronics  were  reviewed  in  Reference 
1  .  The  development  of  tne  device  centered  around  obtaining  a  real  time  elec¬ 
tronic  system  which  would  generate  characteristics  of  any  typical  resistance- time 
curve.  Figure  56  is  a  schematic  of  a  typical  metal  contact  resistance-versu 
time  plot  under  partial  lubrication  conditions  of  sliding  and/or  rolling.  The 
instantaneous  fluctuation  of  the  resistance  is  a  complex  function  of 

(a)  Lubricant  film  thickness,  electrical  properties,  and 

(b)  Contacting  metal  surface  curvatures,  surface  finish  rough¬ 
ness,  and  relative  velocity  between  contacting  surfaces. 

The  Asperitac  provides  a  convenient  t. -'ar.  of  determining  for  a  given  set  of  test 
conditions: 

(a)  The  percent  of  time  spent  in  contact  in  any  lubricated 
metal  interface. 

(b)  The  number  of  contacts  occurring  per  unit  time  across 
the  interface. 

(c)  The  area  of  the  shaded  portion  indicated  in  the  schematic 
of  Figure  56. 

(d)  The  nominal  time-averaged  resistance  present  in  the  con¬ 
tact  zone. 

In  addition,  the  applied  voltage,  center-scale  resistance,  .>nd  resistance 
level  of  detection  can  all  be  preselected  with  the  electronics  used. 

The  percent  of  time  a  lubricated  contact,  spends  in  the  contact  mode  can  be 
obtained  either  digitally  or  in  an  analog  fashion.  The  observed  percent  con¬ 
tact  is  taken  as  the  time  spent  below  the  preselected  detection  rasistance 
relative  to  the  total  time  of  observation.  Figure  57  displays  a  typical  analog 
output  of  percent  contact  on  the  disc  rig  as  a  function  of  both  Hertz  contact 
pressure  and  rolling  speed  when  polyphenyl  ether  is  used  as  a  lubricant.  The 
discs  used  for  the  asperity  work  were  the  same  as  those  used  in  determining 
capacitive  film  thickness  measurements.  The  lower  disc  was  isolated  electri¬ 
cally  from  ground  and  held  at  a  fixed  voltage  level  of  .391  volts  D.C.  for  all 
contact  work  with  the  regulated  power  supply  in  the  Asperitac, 

Since  nominal  film  thicknesses  a  'e  typically  in  the  microinch  region  for  lub¬ 
ricated  contacts,  applied  voltages  in  the  one  volt  range  may  break  down  the 
lubricant  through  electrical  discharge  No  such  breakdown  of  polyphenyl  ether 
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was  observed  up  to  the  1.0  volt  D.C.  available  with  the  power  supply  used. 

Percent  contact  and  asperity  count  rate  was  observed  to  be  constant  for  applied 
voltages  up  to  0.8  volts  across  the  contact  zone. 

The  Asperitac  has  in  its  present  form  a  single  level-crossing  detection  cir¬ 
cuit  and  does  not,  with  a  single  observation,  determine  the  distribution  of 
the  resistance  fluctuation  levels.  This  information  was  determined  by  making 
separate  observations  using  center-scale  resistances  of  10,  100,  1,000,  10,000, 
100,000  ohms.  The  resulting  analog  signals  are  reprinted  in  Figure  58.  Since 
each  detection  level  determines  the  percent  of  time  spent  below  that  resistance, 
a  full  distribution  of  the  fluctuating  resistance  levels  with  time  can  be  deter¬ 
mined  for  a  given  rolling  contact  condition. 

1.  TIME-AVERAGED  RESISTANCE  VERSUS  SPEED 

The  distribution  of  fluctuations  for  any  speed  over  the  test  range  observed 
shows  a  log-normal  dependence  between  resistance  level  and  percent  contc.*-. 

The  cross  plots  for  several  fixed  speeds  taken  from  Figure  58  are  shown  in 
Figure  59.  By  using  the  log-normal  dependence,  one  is  able  to  determine  actual 
time-averaged  resistances  for  every  percent  contact  output  of  the  Asperitac. 

The  percent-contact- time-averaged-resistance  points  are  connected  in  Figure  59 
by  the  dotted  line  shown.  Since  percent  contact  is  observed  over  a  full  range 
from  100  percent  to  zero,  the  functional  dependence  of  nominal  resistance 
averaged  over  time  and  the  full  contact  zone  was  found.  Figure  60  shows  the 
exponential  function  for  a  Hertzian  contact  pressure  of  172,000  psi.  The  ex¬ 
perimentally-observed  exponential  increase  in  resistance  with  speed  turns  out 
to  be  a  sensitive  lubricant  f  .  detector  under  partial  elastohydrodynamic 
conditions . 

At  large  loads  or  low  rolling  speeds,  the  lubricant  film  thickness  will  be  of 
the  same  order  of  surface  roughness  as  the  roiling  surfaces.  Thus,  asperities 
will  have  metallic  interactions  and  part  of  the  applied  load  will  be  shared  by 
interacting  asperities,  Wich  the  description  of  surface  topography  and  of  an 
appropriate  model  for  the  mechanics  of  interactions  between  a  pair  of  asperities, 
such  load  sharing,  real  area  of  contact,  and  electrical  contact  resistance  can 
be  estimated.  A  detailed  analysis  of  the  problem  has  been  described  in  Ref¬ 
erences  1,  16  and  17.  We  shall  present  results  for  typical  ball  bearing  sur¬ 
faces  having  rms  roughness  of  the  order  of  2  microinches.  After  summarizing 
the  surface  topographical  data,  results  of  asperity  load  sharing  will  be  pre¬ 
sented  . 

2 .  SURFACE  TOPOGRAPHY 

Steel  specimens  (AJSI  52100;  in  the  form  of  3  inch  diameter  and  1  inch  thick 
discs  (no  crown)  were  commercially  prepared.  The  rms  surface  roughness  of  less 
than  2  microinches  was  ensured.  Measurements  of  the  surface  topography  con¬ 
sisted  of  obtaining  several  surface  pr  .iles  along  the  axis  of  the  disc  using 
the  Talysurf  Model  4  profilometer  with  a  standard  100  microinch  stylus  (see 
F;gure  61  for  another  set  of  test  surface  tracings).  An  electrical  signal 
proportional  to  the  surface  profile  was  picked  up  from  the  electronic  unit  of 
the  profilometer  and  was  recorded  on  an  FM  tape  recorder.  The  dat  was  then 
digitized  at  a  suitable  sampling  rate  determined  by  highest  frequencies  present 
in  the  data.  A  sampling  interval  of  40  microinches  was  found  to  be  sufficient  for 
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Fig.  59  Percent  Contact  Versus  Contact  Resistance  for  Five  Speed  Leveis 
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Fig.  60  Time-Averaged  Rolling  Contact  Resistance  as  a  Function  of  Speed 
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Nine  Parallel  Talysurf  Tracings  of  Roller  Specimen  Surface 
Finish 


most  surface  profiles  analyzed.  As  will  be  shown  later,  most  surfaces  had  very 
low  frequencies  and  therefore  sampling  at  large  distances  was  necessary  to  ob¬ 
tain  a  reasonable  spectral  distribution  of  the  data. 

Peak  asperity  heights  and  radii  of  curvature  were  estimated  by  using  the  most 
common  three-point  definition  for  a  peak.  Statistical  distribution  of  peak 
heights  is  shown  in  Figure  62.  The  data  is  fitted  by  a  Gaussian  function  and 
the  parameters  of  the  distribution  are  used  in  the  analysis  of  the  surface  in¬ 
teractions  as  described  in  Reference  1.  Radii  of  curvature  at  the  peaks  u 
fitted  by  a  log  normal  function  as  shown  in  Figure  63.  Statistical  correlation 
between  the  heights  and  radii  was  found  to  be  less  than  0.1.  Thus,  tb^  assump¬ 
tion  of  independent  heights  and  curvatures  is  valid. 

Approximation  of  the  three-dimensional  hill  or  summits  by  peaks  given  by  a  sur¬ 
face  profile  has  been  a  very  controversial  question  for  the  last  few  years. 

Nayak  (18)  has  performed  a  detailed  mathematical  analysis  based  on  earlier  work 
(20).  It  is  shown  that  the  differences  between  heights  and  curvatures  of  peaks 
can  be  substantially  different  from  those  of  the  summits.  The  differences  are 
shown  to  be  generally  dependent  on  the  spectral  density  function  of  the  surface 
profile.  If  f(f)  is  defined  as  the  spectral  density  function  where  f  is  the 
frequency,  the  moments  of  ♦(f)  are  defined  as 

to 

u»i  -  J  ♦<£)  f1  df 

o 

It  is  shown  that  an  important  parameter  which  represents  the  differences  be¬ 
tween  profile  and  surface  characteristics  is 

“  =  m0m4/m22 

For  large  a ,  heights  and  curvatures  are  essentially  independent  and  the  differ¬ 
ences  between  profile  and  surface  characteristics  are  expected  to  be  quite  small. 

The  data  of  five  surface  profiles  (Fig.  1)  were  used  to  coiupi te  the  autocorrela¬ 
tion  and  spectral  density  function.  The  recordings  contained  information  di¬ 
rectly  related  to  the  recording  technique  which  was  not  representative  of  the 
local  surface  structure.  The  necessary  components  of  the  data  which  were  re¬ 
moved  before  computing  1r(f)  included  absolute  shifts  in  D.C.  level  signals, 
linearly  increasing  or  decreasing  variations,  and  low  frequency  components  with 
wavelengths  of  the  order  of  the  profile  tracing  length.  The  spectrum  showed  a 
large  peak  in  the  very  low  frequency  region  and  was  close  to  zero  at  higher 
frequencies.  By  sampling  at  large  intervals  the  frequency  band  could  be  re¬ 
duced  to  the  region  of  interest  but  the  upper  cut-off  frequency  introduced  some 
aliasing  error  on  the  trailing  end  of  the  spectrum.  This  is  shown  in  Figure 
64.  The  decay  in  the  autocorrelation  is  also  very  slow;  this  is  shown  separ¬ 
ately  in  Fig.  65,  indicating  dominance  of  low  frequency  waviness  on  the  surface. 

Moments  of  the  spectral  density  functions  were  computed  and  a  was  found  to  be 
close  to  25.  Thus,  the  computed  small  correlations  between  peak  heights  and 
curvatures  is  ptobably  in  agreement  with  Nayak’s  analysis  (18).  Also,  we  are 
justified  in  approximating  characteristics  of  the  surface  by  those  of  the 
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Fig.  o2  Peak  Height  Distribution,  M  =  0.342|i"  and  o  =  1.526u” 
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profile  produced  by  the  talysurf  trace  recordings. 

3.  ASPERITY  LOAD  SHARING 

With  given  statistical  distributions  of  peak  heights  and  curvatures,  the  loads 
supported  by  the  asperities  and  the  size  distribution  of  microcontacts  were 
calculated  using  the  analysis  described  in  Reference  1.  All  asperity  inter¬ 
actions  were  found  to  be  elastic  in  the  range  of  mean  plane  separation  0.5cr.to 
4<T,  where  o  is  the  combined  variance  of  the  peak  height  distribution  of  the 
two  interacting  surfaces.  Identical  surface  topography  was  assumed  on  both 
surfaces  and  elastic  Junction  models  based  on  the  Hertz  solutions  were  used 
to  calc:,  ce  the  real  contact  area  and  load  supported  by  a  pair  of  asperities 
with  given  geometric  interference.  Ouasi- static  sliding  was  allowed  with  a 
constant  interfacial  shear  stress  equal  to  the  flow  stress  («  yield  stress/2) 
of  the  deformating  materials.  The  detailed  analysis  of  the  junction  model  is 
described  in  Reference  17. 

Figure  66  shows  the  loads  supported  by  the  asperities  as  a  function  of  the 
mean  plane  separation  (i.e.,  the  not  .?.)  film  thickness)  for  a  unit  apparent 
area.  The  real  area  of  contact  is  also  plotted.  With  computed  micro  contact 
size  distribution  and  knowledge  of  the  electrical  properties  of  the  materials 
and  expected  lubricant  films,  the  elecgrical  contact  resistance  was  computed. 

For  any  given  geometry  of  interacting  surfaces,  in  our  case  cylindrically 
crowned  discs,  the  Hertzian  contact  area  and  the  expected  lubricant  film  thick¬ 
ness  can  be  estimated.  Dowson  and  Higginson's  formula  was  used  to  compute  the 
film  thickness  as  a  function  of  speed  and  lubricant  properties.  Knowing  the 
film  thickness  and  the  Hertzian  area  (which  will  be  the  total  effective  appar¬ 
ent  contact  area)  the  total  normal  and  friction  load  supported  by  the  inter¬ 
acting  asperities  and  the  real  contact  area  may  be  found  directly  from  Figure 

66.  Solutions  given  in  Figure  66  were  assumed  to  be  valid  for  crowned  discs 
such  as  the  one  used  experimentally  since  all  specimens  were  produced  at  one 
time  in  a  single  lot  with  the  same  techniques.  Thus,  the  loads  shared  by  the 
asperities  may  be  plotted  as  a  function  of  rolling  speed  as  shown  ir  Figure 

67.  The  speed  ranges  for  full  hydrodynamic,  partial  elastohydrodynamic  and 
essentially  dry  contact  are  clearly  evident  from  Figure  67. 

4.  ELECTRICAL  CONTACT  RESISTANCE  AND  FILM  THICKNESS 

For  a  given  film  thickness,  it  is  possible  to  compute  the  size  distribution  of 
microcontacts  and  subsequently  the  electrical  contact  resistance.  Contact  re¬ 
sistance  of  a  single  microccntact  is  given  by  (19) 


where  £  is  the  sum  of  electrical  resistivity  of  the  two  materials,  §  is  the 
tunnel  resistivity  due  to  adsorbed  layers  on  the  mating  surfaces,  and  a^  is 
the  microcontact  radius.  Tunnel  resistivity  5  depends  very  strongly  on  the 
thickness  of  the  adsorbed  layer  (19)  and  may  be  approximately  expressed  as 
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Fig.  66  Average  Normal  and  Friction  Forces  Supported  by  the  Separations 
Between  Mean  Planes  of  the  Interacting  Asperities 
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where 


7.32  x  103  [t  -  (7,2/<J>)  ] 


1.265  x  10 


-  (10/t  er) 


t  *  surface  film  thickness,  (A) 

4>  ■  work  function  for  electron  emission,  (eV) 

ey  »  relative  permittivity  of  the  material  in  the  gap 

The  total  resistance  may  be  obtained  by  summation  over  all  the  microcontacts. 

If  a  model  for  film  thickness  as  a  function  of  speed  is  assumed,  then  it  is 
possible  to  estimate  the  electrical  contact  resistance  if  the  tunnel  resis¬ 
tivity  is  known, 

A  computation  of  tunnel  resistivity  based  on  the  range  of  expected  lubricant 
adsorbed  layers  (as  per  table  in  Reference  20)  showed  that  the  constriction 
resistance  as  compared  to  tunnel  resistance  at  the  microasperities  is  gener¬ 
ally  very  small  and  thus  the  computed  resistance  can  approximately  be  stated 
as  directly  proportional  to  the  tunnel  resistance.  With  f  =  0,5  eV  and 
er  ■  4S  the  values  of  5  as  a  function  of  layer  thickness  t  are  shown  in  Fig¬ 
ure  68, 

For  a  unit  tunnel  resistivity  and  a  unit  apparent  contact  area  the  computed 
resistance  may  be  plotted  as  a  function  of  lubricant  film  thickness  as  in 
Figure  69.  Now  for  any  given  load  knowing  the  Hertzian  or  apparent  area  the 
total  contact  resistance  at  a  given  film  thickness  and  tunnel  resistivity  may 
be  directly  determined  from  Figure  69.  For  the  760  pound,  load,  the  average 
film  thickness  as  a  function  of  speed  was  estimated  using  Grubin's  formula  (10) 
and  thus  the  mean  resistance  was  computed  as  a  function  of  speed  for  three 
values  of  fixed  tunnel  resistivity.  The  results  are  shown  in  Figure  70  along 
with  the  experimentally-determined  exponential  variations  of  resistance  with 
speed , 

The  extremely  high  sensitivity  of  contact  resistance  with  speed,  as  compared 
with  analytical  predictions,  could  be  explained  in  several  ways..  It  is  possi¬ 
ble  that  a  few  scratches  on  the  surface  of  the  discs  may  introduce  highly 
non-Gaussian  characteristics  in  the  surfaces.  However,  it  is  not  clear  as  to 
what  degree  such  deviations  would  change  the  resistance  speed  relationships 
except  that  the  presence  of  a  very  small  number  of  large  asperities  could 
radically  steepen  the  resistance-speed  curves.  Additional  analysis  with  non- 
Gaussian  and  nonisotropic  surfaces  is  needed  to  determine  the  sensitivity  of 
any  such  effects.  Conduction  through  very  thin  absorbed  films  at  interacting 
asperities  is  another  possible  factor  which  could  cause  a  high  sensitivity  of 
contact  resistance  with  speed.  The  resistance  would  be  very  low  when  there 
was  metal  contact,  but  as  the  micro  film  built  up  with  speed  the  interactions 
wcxild  cease  with  a  small  increase  In  speed  and  resistance  would  increase.  When  the  sep¬ 
aration  between  the  mating  asperities  was  in  the  range  of  0  to  0.2  microinch, 
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the  resistances  due  to  a  thin  absorbed  film  would  be  in  the  range  measured. 
Since  this  resistance  Is  a  strong  function  of  the  separations  between  the  as¬ 
perities,  a  high  resistance  speed  sensitivity  would  be  expected. 

The  Asperitac  measurements  may  be  used  to  provide  a  good  indication  of  the 
relative  extent  of  surface  contact,  particularly  in  the  incipit nt  contact  re¬ 
gime.  The  higher  goal,  using  this  device  for  extraction  of  film  thickness  in 
partial  EHD  lubrication,  requires  better  agreement  between  theory  and  experi¬ 
ment  than  has  been  achieved  thus  far.  It  is  anticipated  that  more  analytical 
and  experimental  investigation  along  the  lines  given  above  would  yield  a 
technique  for  determining  lubricant  film  thicknesses  under  partial  elastohy- 
drodynamic  conditions. 
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SECTION  VII 
CONCLUSIONS 


The  traction  data  obtained  with  the  MIL-L-7808  oil  showed  both  quantitative  and 
qualitative  dissimilarities  with  the  5P4E  data,  even  under  conditions  where 
conventional  viscosities  and  pressure  coefficients  of  viscosity  were  similar . 
Order  of  magnitude  differences  in  tractions  measured  under  these  conditions 
show  that  great  caution  must  be  taken  when  specifying  lubricant  properties  in 
that  conventional  properties  do  not  appear  to  determine  lubricant  traction  be¬ 
havior.  In  general,  the  MIL-L-7808  oil  produced  markedly  lower  tractions  at 
all  temperatures  studied  and  showed  a  significant  temperature  dependence  which 
was  not  exhibited  by  the  polyphenyl  ether.  A  hypothetical  viscosity  model  has 
been  put  forth  (as  was  done  for  polyphenyl  ether)  which  mathematically  repre¬ 
sents  the  traction  data  over  the  ranges  of  loads,  speeds,  and  temperatures 
studied  experimentally.  It  was  necessary,  however,  to  allow  for  time  effects 
when  obtaining  this  correlation  since  the  variations  in  traction  with  rolling 
speed  are  much  larger  than  can  be  accounted  for  by  predicted  variations  in  film 
thickness  alone. 

Film  thicknesses  were  measured  by  means  of  optical  and  capacitance  techniques. 
The  optical  measurements  provided  the  high  degree  of  accuracy  observed  by  other 
investigators  who  have  used  this  technique.  In  particular,  it  showed  the  trail¬ 
ing  edge  dip  which  for  the  conditions  of  highly-elliptical  contact  extended 
across  the  back  but  did  net  appear  to  come  around  the  sides.  The  optical  mea- 
suremenets  obtained  with  one  cylindrical  transparent  quartz  disc  and  one  pol¬ 
ished  steel  disc  at  fairly  low  pressures  and  speeds  seem  to  be  in  reasonably 
good  agreement  with  Dowson  and  Higginson’s  isothermal  theory.  This  agreement 
has  been  borne  out  by  optical  measurements  obtained  by  Cameron,  particularly 
regarding  the  insensitivity  of  film  thickness  to  load. 

Capacitance  measurements,  on  the  other  hand,  indicated  similar  speed  variations 
as  predicted  by  elastohydrodynamic  theory;  but  showed  considerably  higher  var¬ 
iations  with  load  at  the  relatively  higher  loads  and  speeds  that  were  studied. 
These  load  variations  are  borne  out  in  part  by  the  Battelle  X-ray  data  where 
even  stronger  variations  in  load  were  indicated  under  these  conditions.  The 
full  cause  of  these  variations  is  still  at  this  time  not  fully  understood, 
although  it  is  believed  that  the  Battelle  measurements  have  been  at  least 
qualitatively  substantiated. 

Asperitac  measurements  indicated  partial  elastohydrodynamic  contact  to  start 
at  speeds  somewhat  below  100  inches  per  second  which  corresponds  to  predicted 
film  thicknesses  below  20  microinches.  More  work  is  needed  in  this  important 
area  to  fully  correlate  Asperitac  data  with  asperity  interaction  analyses  so 
that  effective  film  thickness  data  in  the  partial  elastohydrodynamic  regime 
can  be  established  and  used  to  predict  partial  elastohydrodynamic  performance 
characteristics.  It  is  believed  that  this  overall  approach,  if  continued, 
would  be  quite  fruitful. 
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APPENDIX  I 


PHYSICAL  PROPERTIES  OF  MIL-L-7808 
LUBRICANT  USED  IN  THIS  INVESTIGATION 


Index  of  Refraction  -  1.395 

Density  -  See  Figure  71 

Viscosity  -  See  Figures  16,  72,  and  73 

Bulk  Modulus  -  See  Figure  74 

Pressure-Viscosity  Coefficient  -  See  Formula  in  Text  (Page  22) 
Temperature-Viscosity  Coefficient  -  See  Formula  in  Text  (Page  22) 
Thermal  Conductivity  -  .056  Btu/hr/ft/°F 
Specific  Heat  -  .5  Btu/lb-°F 
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Fig.  71  Density  as  a  Function  of  Temperature  and  Pressure  for  0-67-29 
(MIL- L- 7808) 
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Kinematic  Viscosity  as  a  Unction  of  Temperature  and  Pressure  for 
0-67-20  (MIL- L- 7808) 


APPENDIX  II 


TRACTION  DATA  SUMMARY  FOR  TURBINE  ENGINE  LUBRICANT 
MIL-L-7808G  AIR  FORCE  CODE  0-67-23  (SEE  APPENDIX  I  FOR  IDENTIFICATION) 


Test  Condition  Summary 


MTI  Computer  Data 

Rolling  Speed 

Hertzian  Contact 

Lubricant  Test 

Bank  Curve  Number 

in/ sec 

Pressure,  psi 

Temperature, 

62 

448. 

100,000. 

78. 

63 

448. 

120,000. 

78. 

64 

448. 

130,000. 

74. 

65 

448. 

130,000. 

78. 

66 

448. 

100,000. 

100. 

67 

448. 

120,000. 

100. 

68 

448. 

130,000. 

100. 

69 

448. 

130,000. 

100. 

70 

448. 

100,000. 

120. 

71 

448. 

120,000. 

120. 

72 

448  o 

130,000. 

120. 

73 

448. 

130,000. 

120. 

74 

448, 

140,000. 

120. 

75 

448. 

100,000. 

140. 

76 

448. 

120,000. 

140. 

77 

448. 

130,000. 

140. 

78 

448. 

130,000. 

140. 

79 

448. 

140,000. 

140. 

80 

895. 

100,000. 

78. 

81 

895. 

120,000. 

78. 

82 

895. 

130,000. 

78. 

83 

895. 

140,000. 

78. 

84 

895. 

100,000. 

101. 

85 

895. 

120,000. 

101. 

86 

895. 

130,000. 

101. 

87 

895. 

140,000. 

101. 

88 

895. 

100,000. 

120. 

89 

895. 

120,000. 

120. 

90 

895. 

130,000. 

120. 

91 

895. 

140,000. 

120. 

92 

895. 

100,000. 

140. 

93 

895, 

120,000. 

140. 

94 

895. 

130,000. 

140. 

95 

895. 

140,000. 

140. 

96 

1351. 

100,000. 

75. 
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Test  Condition  Summary  (Continued) 


MTX  Computer  Data 

Rolling  Speed 

Hertzian  Contact 

Lubricant  Test 

Bank  Curve  Number 

in/ sec 

Pressure,  psi 

Temperature,  •' 

97 

1351. 

120,000. 

75. 

98 

1351. 

130,000. 

75. 

99 

1351. 

140,000. 

75. 

100 

1351. 

100,000. 

102. 

101 

1351. 

120,000. 

102. 

102 

1351. 

130,000. 

102. 

103 

1351. 

140,000. 

102. 

104 

1351. 

100,000. 

120. 

105 

1351. 

120,000. 

120. 

106 

1351. 

130,000. 

120. 

107 

1351. 

140,000. 

120. 

108 

1351. 

100,000. 

140. 

109 

1351. 

120,000. 

140. 

110 

1351. 

130,000. 

140. 

111 

1351. 

140,000. 

140. 

112 

1351. 

100,000. 

150. 

113 

1351. 

120,000. 

150. 

114 

1351. 

130,000. 

150. 

115 

1351. 

140,000. 

150. 

116 

1351. 

100,000. 

160. 

117 

1351. 

120,000. 

160. 

118 

1351. 

130,000. 

160. 

119 

1351. 

140,000. 

160. 

120 

1351. 

150,000. 

160. 

121 

1822. 

100,000. 

80. 

122 

1822. 

120,000. 

80. 

123 

1822. 

130,000. 

80. 

124 

1822. 

140,000. 

80. 

125 

1822. 

150,000. 

80. 

126 

1822. 

100,000. 

100. 

127 

1822. 

120,000. 

100. 

128 

1822. 

130,000. 

100. 

129 

1822. 

140,000. 

100. 

130 

1822. 

150,000. 

100. 

131 

1822. 

100,000. 

140. 

132 

1822. 

120,000. 

140. 

133 

1822. 

130,000. 

140. 

134 

1822. 

140,000. 

140. 

135 

1822. 

150,000. 

140. 

NOTE:  Test  Disc  Diameter:  3.00  in. 

Surface  Finish  :  ~  2  microinch  rms 

Crown  Radius  :  36  inches 

Material  :  AISZ  52100  bearing  steel 

Owes  1-61  :  Contain  polyphenyl  ether  data  (Ref.  1) 
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Fig.  77  Computer  Plotted  Experimental  Traction  Versus  Slip  Data 
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Fig.  82  Computer  Plotted  Experimental  Traction  Versus  Slip  Data 
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Computer  Plotted  Experimental  Traction  Versus  Slip  Data 
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Fig.  87  Computer  Plotted  Experimental  Traction  Versus  Slip  Data 
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APPENDIX  III 

ELASTIC  DISTORTION  ANALYSES  USED  FOR 
STATIC  OPTICAL  CALIBRATION  AND  CAPACITANCE  CALIBRATIONS 


The  ability  to  determine  the  shape  of  the  full  exterior  Hertz  contact  profile 
in  elliptic  contact  is  very  useful  in  calibrating  various  electrical  and  op¬ 
tical  film  thickness  in  concentrated  contact  elements.  The  basic  equation 
for  determining  the  elastic  distortion  profiles  were  first  put  forth  by  Hertz. 
These  equations  were  presented  in  convenient  computational  form  in  Reference 
(l)and  were  used  to  obtain  comparisons  between  theoretical  and  measured  opti¬ 
cal  interference  fringe  patterns  for  the  particular  disc  geometry  used  in  the 
investigation.  This  appendix  will  be  devoted  to  the  presentation  of 
dimensionless  curves  and  methods  for  extending  these  predictions  to  any  ellip¬ 
tical  contact  geometry  and  to  presenting  the  methods  used  in  obtaining 
capacitance  film  thickness  calibrations. 


Formulas  for  Exterior  Hertz  Profile  in  Elliptic  Contact 

The  formulas  presented  in  Reference  1  are  based  upon  those  originally  put 
forth  by  Hertz  and  are  repeated  here  for  completeness. 


A  plan  and  elevation  view  of  an  elliptical  Hertz  contact  is  shown  in  Figure 
III-l.  The  exterior  Hertz  profile  h  is  expressed  as  a  function  of  the  cr or¬ 
dinate  x  and  y  in  terms  of  the  elastic  modulus  and  Poisson's  ratio  of  contact¬ 
ing  bodies  1  and  2  (denoted  E^,  v^,  and  v  respectively),  the  applied  load 

W  and  the  contact  ellipse  dimensions  a  and  b  in  terms  of  the  dimensionless 
quantities  <i>,  5,  and  n: 


<5> 


2E*a 

3W 


h,  5 


x 

b  » 


8  = 


b 

a 


where 


E1 


(1  -  v^)  (1  -  v22) 


nE, 


ttE, 


The  Dimensionless  separation  $(5,0)  is  determined  for  Hertz  contact  theory 
which  may  be  put  in  the  form 

♦<C.n>  *  J252  +  J3n2  -  J1  (A3.1) 

where 


ru 


j 


dC 


(1  + 


;V/2 


<1 


sW71 


(A3. 2) 
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and  u  Is  determined  from  the  relationship 


_2  2 

5  ,  n 

''  "15  + - 121? 

1  +  u  l  +  H  u 


(A3. 5) 


The  integrals  J, ,  J9,  and  Jo  may  be  expressed  in  terms  of  the  elliptic  inte- 
'  —  •*  "3(0,  ■  ‘  J  '  “ 


grals  F(0,  k)  and  E(0,  k)  as  follows: 


F(0,  k). 


J2  " 


1Zr  [E(0,  k)  -  B2F(0,  k)] 


a  -  r> 


and 


(1 


■hr  IF«- k>  * E«> k»  +  n  ,  Ju,i/2 ,,  ,  -KUi 

-  3  )  (1  +  M  )  (1  +  u  .) 


where 


'  -2  -1  - 
y  1  -  B  and  0  »  tan  u 


and  F(0,  k)  and  E(0,  k)  denote  the  incomplete  elliptic  integrals  of  the  first 
and  second  kind  defined  as 

A  f  A  . _ 

F(0,  k)  «  /  -  -  =P==r  ,  and  E(0,  k)  -  /  VI  -  k2sin20*  d0' 

V 1  -  k2sin20!  o-' 

The  aspect  ratio  of  the  contact  ellipse  3  is  related  to  the  relative  curvature 
parameter  R/R’  in  accordance  with  the  equation 


2  [K(k)  -  E(k) ]  ,  _R 

E(k)  -  82K(k)  R' 


(A3. 6) 
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The  dimensional  Hertz  contact  parameters  are  calculated  from  the  relationship 

(A3. 7) 


rl  IH  _JR1a1/3  " 
*2  E’  R  +  R* '  n 


and 


b/8 


P  .  1 JL 

HZ  2  nab 


(A3. 8) 


(A3. 9) 


The  relationships  between  8,  n,  and  R/R'  appear  frequently  in  the  literature 
but  are  repeated  here  in  Figure  III-2  for  completeness. 

The  dimensionless  exterior  profile  for  elliptic  contact  as  obtained  from  Equa¬ 
tions  (A3.1)  through  (A3. 5)  does  not  appear  in  the  literature  and  results  of 
computation  for  computing  exterior  Hertz  separation  profiles  are  given  here 
in  the  form  of  values  of  4>  on  the  x  and  y  axes  respectively  in  Figures  III-3 
and  III-4. 

The  family  of  curves  in  Figure  III-3  is  of  the  form  versus  4>(£o,  0)  for 
various  values  of  R/R* ,  where  denotes  values  of  £  when  n  *  0. 

Similarly,  Figure  III-4  shows  values  of  versus  $(0,  nD)  where  nQ  corresponds 
to  values  of  n  when  5=0.  For  a  given  dimensionless  separation,  ?,  one  may 
calculate  the  values  of  £0  and  ri^  corresponding  to  the  extreme  points  of  the 
closed  curve  representing  the  locus  points  of  constant  separation  4>  (we  will 
refer  to  these  curves  subsequently  as  fringe  lines). 

Complete  fringe  lines  for  values  of  n0  13  1,  2,  3,  4,  and  5  and  R/R'  *  1/24 
(corresponding  to  3  inch  diameter  discs  with  36  inch  crown  radii)  are  shown  in 
Figure  I1I-5.  These  are  plots  of  n  versus  C  at  constant  £.  The  curves  are 
very  nearly  elliptical  in  nature  and  in  the  extreme  cases  n0  ■  1  and  hQ  -*•  ", 
the  curves  approach  those  of  true  ellipses.  These  represent  the  Hertz  contact 
ellipse  and  rigid  body  separation  fringe  lines  respectively.  Good  approxima¬ 
tions  to  the  fringe  lines  can  be  obtained  with  the  ellipses  having  major  and 
minor  semi-axes  of  nQ  and  £0.  Hence,  lines  of  constant  separation  $  may  be 
approximated  by 

(^-)2  +  <7p)2  -  1  (A3. 10) 

S  o 

where  and  may  be  obtained  as  functions  of  ♦  from  Figures  III-3  and  III-4. 
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Fig.  93  Hertz  Contact  Parameters 
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Fig.  96  Dimensionless  Fringe  Lines  for  R/R'  =  1/24 


The  maximum  error  for  the  test  case  under  consideration  (R/R*  ■  1/24)  is  about 
6  percent  and  occurs  when  n0  *  1.2.  The  elliptical  and  true  fringe  lines  for 
this  case  are  shown  in  Figure  24  appearing  in  Section  III  of  this  report. 

This  approximation  will  be  used  in  the  following  section  when  calculating 
capacitance  and  will  provide  a  considerable  savings  in  computer  costs  with 
little  error  in  accuracy. 

Capacitance  Calculations 

In  order  to  calculate  capacitance  it  is  necessary  to  obtain  an  effective  area- 
to-displacement  relationship  (A'/h) .  If  h  is  to  be  assumed  constant  on  ellip¬ 
tical  fringe  lines,  then  (A/h)  may  be  calculated  by  dividing  differential 
elliptical  bands  by  the  separation  at  the  band  and  integrating  over  the  entire 
space  bounded  by  the  discs.  The  elliptical  bands  under  consideration  are 
shown  schematically  in  Figure  III-6.  They  fall  into  two  categories  correspond¬ 
ing  to  whether  the  fringe  line  falls  completely  between  the  teo  discs  as  shown 
in  Figure  Ill-6a  or  is  truncated  by  the  end  of  the  discs  having  axial  half 
length,  L,  as  shown  in  Figure  III-6b. 

The  area,  A,  bounded  by  a  line  of  constant  separation  approximated  by  the 
ellipse  having  semi-axes  xq,  yQ  is  given  by 

A  »  J  x  y  ,  y  <  L  (A3. 11) 

o  o  o 


When  these  ellipses  go  beyond  the  edge  of  the  disc  having  half  thickness,  L, 
then  the  area  A  bounded  by  the  solid  lines  in  Figure  III-6b  is  given  by 

A  ■  2  “o  yo  ‘M/1  -  (r)2  +  (7-)].  y0  >  1-  (A3 . 12) 

For  any  given  value  of  R/R*  and  $  one  can  determine  x  “  b 5  within  the  use 
of  Figure  III-3  and  yQ  *  anQ  within  the  use  of  Figure°III-4,  and  the  necessary 
crossplots  may  be  obtained  to  provide  x0  as  a  function  of  yQ  which  are  presen¬ 
ted  in  the  dimensionless  form  of  nQ  versus  ?0  in  Figure  III-7, 


If  we  adopt  the  Grubin  assumption  that  the  shape  of  the  film  thickness  profile 
follows  that  of  the  exterior  Hertz  profile,  then  the  EHD  film  thickness  profile 
h  may  be  expressed  mathematically  as 


h 


hl  + 


3W  t 


2E’a 


(0,  yQ/a). 


y  <  a 

o 

y  >  a 

o 


W  +  0 


(A3. 13) 


together  with  the  constraint  that  h  is  assumed  constant  on  the  fringe  lines 
given  by 

(^*)2  +  (~ )2  -  1  (A3. 14) 

yo  o 
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The  indicated  area  derivative  may  be  obtained  by  differentiation  of  Equations 
(A3. 11)  and  (A3. 12)  and  is  given  as 


ir(x  +  y  -~-) 
o  7o  dy 

o 


2L(-~  -  -2-) 
d*0 


,  yD  <  l 


T  2  dx  1  T 

<“>  +2ydrsin 

O  0  o 


(A3. 17) 


where  dx  /dy  -  6d  /d?  may  be  obtained  by  differentiation  of  Figure  III-7. 
0  0  o  o 


When  the  discs  are  unloaded.  Equations  (A3. 14)  through  (A3. 17)  remain  the  same; 
but  the  elastic  distortion  profile  is  replaced  by  the  rigid  body  profile  which 
may  be  expressed  by  mathematically  replacing  Equation  (A3. 13)  with 


,  W  -  0 


(A3. 18) 


where  is  related  to  yQ  by 

xq  =  V R/R'  yQ  ,  W  *  0  (A3. 19) 

It  should  be  noted  that  a  smooth  transition  does  occur  between  lightly  loaded 
and  unloaded  conditions. 

Equation  (A3. 16)  may  be  readily  integrated  numerically  by  Gaussian  Quadrature 
to  obtain  capacitance  computation  in  the  form  of  (A/h)  as  functions  of  h^  and 
W.  The  results  of  these  computations  are  presented  and  described  in  the  main 
text  of  this  report. 
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APPENDIX  IV 

DISTORTION  FACTOR  DERIVATION  FOR  INTERFERENCE 
FRINGE  SPACING  IN  ROLLING  DIRECTION  OF  CONTACT  ELLIPSE 


The  geometric  arrangement  depicted  in  Figure  IV-1  shows  an  exaggerated  example 
of  how  i  fringe  image  observed  to  be  at  S'  is  offset  from  the  actual  position 
S  in  the  contact  area.  The  dijc  is  cut  from  a  cylinder  with  a  center  bore 
concentric  to  the  outside  diameter  of  the  disc.,  Both  inner  and  outer  surfaces 
have  been  optically  polished.  The  outside  diameter  has  been  chromium  coated 
and  provides  a  semi-transparent  surface  which  produces  interference  patterns 
of  the  contact  region  gap. 

Let  S  define  the  actual  distance  to  any  fringe  in  the  rolling  direction  along 
a  perpendicular  line  from  the  majo'-  axis  of  the  contact  ellipse.  Let  S’  be 
the  observed  distance  from  the  image  plane  to  any  fringe  along  the  same  per¬ 
pendicular.  For  a  lixed  disc  geometry  and  index  of  refraction,  the  distortion 
factor,  f,  between  observed  images  and  actual  is 

fjj  -  —  f  (" i  R^»  nQ»  L^) 

The  microscope  arrangement  is  such  that  L-^  >>R->  and  for  all  fringing  makes 
0  <  2  The  contact  illuminating  light  can,  therefore,  be  assumed  as  parallel 
to  the  axis  line  cf  the  test  disc 

Since  the  microscope  is  in  air,  Snell’s  Law  states 

sin  0  =  n^,  sin  V  (A4,l) 


Small  angle  approximations  provide  that 

AS 


sin  y  =  tan  y  = 


Rrr 


sin  G=0= 


S’ 


Also,  from  the  diagram  we  note 


S+AS 

k. 


JLL 

r 


Substituting  Equations  (A4..2)  and  (A4„3)  into  (A4d)  we  get 

c»  Vr 

AS  =  (--)  (-r-) 

r  nG 


(A4„2) 

(A4-.3) 


(A',.4) 


(A4 -5) 
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-Axis  Rolling  Direction 


Fig.  99  Optical  Distortion  Off-Set  Diagran 


Evaluation  of  Equation  (A4.4)  with  this  value  of  AS  results  in 


For  the  test  setup  evaluation  yields 

«  1.5  inches 

r  =  1.00  inches 

rG  -  1.458 

«  .8641 


(A4.6) 


In  other  words,  the  visually-observed  fringe  location  will  be  .8641  times  the 
actual  location  from  the  major  axis  line  of  the.  elliptical  co-.tact  pattern. 
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f  APPENDIX  V 

i 

n* 

t  EQUIPMENT  LIST 

i' 

•  Capacitance  Film  Thickness  Measurement 

'  •  MTI  Rolling  Disc  Machine  (described  in  detail  in  Reference  1) 

•  Wayne-Kerr  autobalance  precision  bridge  Model  B331  MK11 

•  Two  mercury-wetted  rotating  contacts 

•  .ae  Mechanical  Technology  Incorporated  differential  tachometers 
(one  side  only  for  absolute  speed  output);  Drawing  No.  10D000068 

•  One  Hewlett-Packard  x-y  recorder  Model  70Q4A 

•  One  Mitutoyo  calibration  tester  micrometer  Model  (521-102) 

Optical  Film  Thickness  Meam.  err.ent 

•  One  Nikon  microscope  with  IX  objective  and  15X  eyepiece  lenses 
and  colored  filters 

•  Two  Dell  Optics  cylindrical  chromium-coated  quartz  discs  per  MTI 
Sketch  No.  SK-A-4211 

•  One  Pentax  35mm  camera  with  microscope  mount 

•  One  three-way  directional  microscope  holder 

•  One  45-degree  extension  tube  mirror  holder 

•  One  Honeywell  65C  stroboscope 

•  Fo tonic  sensor 

•  Digital  Monsanto  counter 

X-Ray  Film  Thickness  Measurement 

•  Harshaw  Chemical  Company  amplifier  and  pulse  height  discriminator 
Model  N-301 

•  Harshaw  Chemical  Company  regulated  power  supply  Model  N-^Ol 

•  Harshaw  Chemical  Company  linear  count  rate  meter  Model  N-730 

109 

•  One  New  England  Nuclear  Cd  source  Serial  No.  D-lll 

109 

•  One  Nuclear  Associated  Cd  calibration  source  Cac.  No.  67-405A 

•  One  Baird-Atomic  X-ray  safety  survey  meter  Model  420 

•  Bendix  safety  pocket  dosimeters  Model  1200-MR 

•  One  lead  container  MTI  design  Drawing  No.  245-D-035 
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